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ABSTRACT 

This thesis describes the design and construction of a wide band 
(60-3060 Hz), three phase (including zero sequence), non-transposed, 
lumped element, power transmission line model. The model consists of 
from 15 to 65 "three phase L sections" - depending upon the length of 
the actual line being modeled - and includes trimming features to model 
both horizontal and vertical configurations of lines in the 345 kV to 
1500 kV range. 

To predict the performance of the model and determine the number 
of sections required to model a given line, the poles of the driving 
point impedance of various numbers of L-sections in series is determined 
by computer and compared with the corresponding poles of the actual line 
being modeled. For either an open or short circuited termination, the 
performance specification is that in the range from 60 to 400 miles long, 
the percent error between the poles of the model and those of the actual 
line shall not exceed 5 percent. This criterion results in the require- 
ment for 15 to 65 L-sections. Thirty sections are of 4 mile "length" 
and thirty five are of 8 mile "length". 

With the number and "length" of the model sections determined, 
electrical parameter data for a number of typical extra high voltage 
(EHV) lines are suitably scaled to select a base value of model inductance 
and capacitance for the 4 mile and 8 mile sections. The effects of non- 
transposition are modeled by incorporating two mutual inductors between 
the phases. Zero-sequence inductance is modeled by a specially designed 
low Q inductor whose inductance decreases appropriately with increasing 
frequency. 

To ensure that the overall quality factor of the model remains 
greater than or equal to that of actual lines for all frequencies between 
60 Hz and 3060 Hz; low loss, powdered molybdenum-nickel-iron toroid cores 
are selected for the inductors, and low dissipation factor polystyrene 
capacitors form the basic capacitance of each section. With the designed 
quality factor of the model exceeding that of most EHV lines in operation, 
the Q of the model is brought into correspondence with actual lines (at 
60 Hz) by suitably inserting resistance into each section. 
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CHAPTER 1 
INTRODUCTION 

U Transmission System Simulators 

The three phase power transmission line model to be designed 
will form part of a power transmission system simulator (TSS) for use 
in the M.I.T. Electric Power Systems Engineering Laboratory. Knowledge 
of the transient behavior of power system networks is essential to the 
system designer and planner in order to correctly dimension or rate 
the various components of a power system. Circuit breakers and 
insulators, for example, cannot be correctly sized without some 
knowledge of the magnitude and duration of system fault currents or 
switching overvoltages. With 765 kV systems already in operation, 
and 1100 kV and 1500 kV systems proposed for the future, it is 
particularly important to know the magnitude and shape of over- 
voltage transients caused by routine switching operations. Models 
for 1100-1500 kV lines still in the design stage are of obvious value, 
and the very limited access to extra high voltage (EHV) lines in 
operation makes a model of these lines highly desirable. 

1.2 General Model Parameters and Constraints 

A transmission line model (t.l.m.) is only one part of a TSS. 
Most simulators include, in addition, models of power system generators, 
transformers, compensating reactors, lightning arresters, circuit 
breakers, ground fault neutralizers, etc. In the present t.l.m. 
design, the requirement for interface with other TSS components places 
certain constraints upon the model. In particular, compatibility 
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with existing TSS transfonners dictates that the nominal operating 

voltage of the t.l.m. be 138V line-to-line and that nominal phase 

current be 38.6 ma. Other required features for the t.l.m. are that 

it be capable of modeling vertical, horizontal, transposed, and non- 

trans posed EHV line configurations from 60 to 400 miles in length, 

have quality factors exceeding those of existing EHV lines from 60 

to 3060 Hz, model the variation of zero sequence inductance with 

frequency over the range 60 to 3060 Hz, and model lines in the 

voltage range from 345 kV to 1500 kV. In addition, the positive 

sequence inductance and capacitance of the t.l.m. must remain 

essentially constant over the range 60 Hz to 3060 Hz, and inductance 

must remain constant with varying current in the range from 0 to at 

least one ampere-thus ensuring constant model inductance up to 

25 p.u. fault current. 

1.3 Features and Scope of Design 

One of the more unique features of this t.l.m. design is 

that it is to be capable of modeling non-trans posed power 

transmission lines. Numerous transmission line models in simulators 

1 -4 

are presently in operation , but all in the U.S., at least, are 
transposed models. A literature search covering the period subsequent 
to 1960 disclosed only a limited literature on t.l.m.'s - all 
relating to transposed models. Many actual lines are, in fact, not 
transposed. Previously, with relatively lower voltage transmission 
lines, it was assumed that the small effects of non transposition 
could be ignored. Non-transposed lines were mathematically and 
physically modeled as though they were transposed, and model results 
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v/ere considered to agree satisfactorily with actual lines. Recent 

5 

studies on EHV lines indicate, however, that the effects of non- 
transposition can no longer be ignored. This thesis is an attempt 
to design a t.l.m. which will include the effects of non-transposition 
and document the design in sufficient detail that persons other than 
the designer can construct and operate the model. 

1.4 Extensions of the Design Concept 

The t.l.m. design discussed in this thesis is intended 
specifically for the modeling of electric utility power transmission 
systems for which the ratio (per unit length) at 60 Hz is of 

the order of 230,000-270,000 rti^ (X^ = positive-sequence capacitive 
reactance; X^^ = positive-sequence inductive reactance). The design 
concept, however, is not at all limited to utility lines. As will 
become clear in later sections of the discussion, any band of 
Xc/X^ ratios could have been provided - had it been desirable or 
necessary. With the proper selection of model components, any audio 
range power transmission system involving inductive and capacitive 
coupling could be accommodated. Thus, the model concept is equally 
suited to modeling, for example, under ground power cables, long 
ship-board or industrial plant cable runs or long synchro cable runs, 
etc.. In addition, the mathematical approach developed in Chapter 3 
for the solution of the impedance polynomial has application in any 
area where the impedance concept has validity - including, for 
example, the mechanical impedance in mechanical vibrations problems. 
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CHAPTER 2 

MODEL CONFIGURATIONS 



2.1 General Non-Transposed Models 

A section of a general, lumped-element, three-phase 
transmission line is shown in Figure 2.1. The electrical para- 
meters of tne model follow directly from the inductance, 
capacitance, and resistance matrices of the actual line. With the 
Ilia trices of the line defined as 



and 



[L] = 



[C] = 



[R] = 



l-aa 


'-ab 


Lac 


Lab 


^bb 


Lbc 


Lac 


4c 


Lee 


Laa 


-4b 


-Lac 


-Lab 


4b 


-Lbc 


-Lac 


-4c 


Lee 


'*aa 


4b 


\c 


"ab 


4b 


"be 


"ac 


4c 


"cc 



(2-1) 



tne correspondence between line and model inductance and capacitance 
parameters is as follovys (See Appendix A.l): 



'-aa ^11 ^ '-00 



■bb 



“cc 



^22 ^00 
^33 ^00 



^ab " ^12 
^ac " ^13 



•-ab " ^12 ^00 

^bc ' 43 4o 
^ac " ^13 4o 



^aa " *'12 ^13 '^lO 

^bb *'12 43 ^ 4o 



( 2 - 2 ) 



4c " ^23 



^cc " ^13 ^23 4o 
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The resistance Rqq in the model represents only an approximate model- 
ing of line inter-phase resistances. Rather than providing three 
separate inter-phase resistance circuits in the model, each of the 
three line inter-phase resistances is approximated by the average 
of the three i .e. , 



''' 'V .. ^ab ^ ^bc ^ ^ac _ 



^ab ~ ^bc ” ^ac 



= Rr. 



(2-3) 



So that. 



% 

[R] = 



^aa ^m 

*^bb ^m 
^in ^m ^cc 



(2-4) 



Comparing line equations vntn those of Appendix (A.l), it is found 
that: 



•^00 = Rm 

*^aa " ^11 *^m ^2. 

*^bb 

^cc " *^33 '^1 

The line parameters given in equation (2-1) include the effect 
of ground wires, so that the model parameters also include ground 
wire effects. The mutual inductances shown in Figure 2.1 represent 
"ideal" mutual inductors adding no self inductance to the lines. 

The value of the mutual inductances is nominally about ten percent 
or less of the self inductances. 

The circuit of Figure 2.1 can be transformed by simple 
algebra (See Appendix A. 2) to the configuration of Figure 2.2. 

The circuit of Figure 2.2 is as general as that of Figure 2.1 but 
has one less mutual inductor. In Figure 2.2 the transformed 
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quantities, in terms of Figure 2.1, are: 
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(A-4) 



If the line being modeled is of horizontal configuration; then for 
both circuits: 
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The model section shown in Figure 2.2 (hereafter referred to as 
"the model") is the configuration selected for the design. While 
it is by no means the only configuration possible, this model has 
a number of advantages over other possible models of equal generality. 
First, it should be pointed out that in order to economize on the 
number of cores to be purchased, the mutual inductors are wound over 
tne self inductors in the line. The configuration of Figure 2.2 
requires four cores, and incorporates the feature that no self- 
inductor will be required to physically accept more than one mutual 
winding (the self winding leaves sufficient window area for only one 
mutual). Second, the configuration in Figure 2.2 ensures that no 
phase line is required to incorporate the resistance of more than 
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one mutual winding. (The resistance of each phase line must be 
kept to a minimum to keep Q high.) Third, when modeling horizontal 
lines, the mutual inductances of the Figure 2.2 configuration are 
equal -thus incorporating into the model the symmetry v;hich actually 
exists in the line being modeled. (This would not be true of some 
other non-transposed models.^) 

2.2 Transposed Models 

As mentioned in Chapter 1, it is common practice to simplify 
the analysis of non-transposed lines by assuming that they are 
transposed. In most non-transposed lines the values of the three self 
inductances of the phase conductors are approximately the same, and the 
values of all itJter-pnase parameters are not too different from one 
another. As shown in Appendix A. 3, the transformation to a transposed 
view-point "averages" the line parameters and "re-distributes" them 
uniformly among the three phases. The result is a balanced, symmetrical 
model which for many applications is a good approximation to what is 
in reality a non-transposed line. A schematic of a transposed line 
model section resulting from the transformation of Appendix A. 3 is 
shown in Figure 2.3. In terms of the general model of Figure 2.1, 



the parameters of Figure 2.3 are: 
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^0 
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2.2.1 Symmetrical Component Equivalent Circuits for Transposed Hodel 
In Appendix A. 3 it is shovjn that the transposed line model 
can be reduced to two symmetrical component models. The results are 
repeated here for continuity. With the definitions 



^P.S. 

^z.s. = 



‘^Z.S. ■ ^^'00 

'-p.s. " ‘-S ■ 41 
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the transposed line equations (A-7) - with positive sequence currents 
only- reduce to: 
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And The Current Law equations (A-11) reduce to: 
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These results show that for positive (or negative) sequence voltages 
and currents, the transposed line may be viewed as three identical 
single phase circuits (the 120° phase difference understood). 

Similar conclusions apply to zero sequence voltages and currents 
where it can be shown (Appendix A. 3) that under these conditions 
the transposed line equations reduce to: 
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and 
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Thus, the non-transposed model can bd approximated by a transposed 
model, and the transposed model can, in turn, be analyzed into two 
symmetrical component equivalent circuits - one for positive-sequence 
and another for zero-sequence. 
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CHAPTER 3 

NATURAL FREQUENCY ANALYSIS 

3.1 Methods and Criteria for Predicting Model and Actual 
Line Correspondence 

As discussed in Chapter 2 and Appendix A, the transposed 
line can be analyzed into two, symmetrical component, single-phase 
equivalent circuits. A single-phase L-section model as "seen" by 
positive (or negative) sequence voltages and currents is shown in 
Figure 3.1a. A similar L-section model as seen by zero-sequence 
quantities is shown in Figure 3.1b. Figures 3.1a and 3.1b reflect 
the basic "L" configuration iriiplied by the three-phase models of 
Figures 2.1, 2.2 and 2.3, but "ir" or "T" section models are equally 
valid. In fact, under certain circumstances and with certain loads, 
the "tt" or "T" configurations may have advantages over the "L" 
models - these circumstances will be discussed more fully in later 
sections. 

One of the purposes for which the transposed line 
approximation can be used is that of establishing a "math-model- 
predictor" for co(nparison of an actual line and its proposed models. 
Some of the early questions in the design of a lumped element 
t.l.m. are; How many sections are required? How "accurate" is the 
model with ten-twenty-forty sections? Should the sections be "L" 
or "tt"? What should be the "length" of each model section? For 
what frequency band must the model "correspond" to the actual line? 
Exactly what is it that should "correspond" between model and line- 
to what "accuracy"? A method of theoretically comparing an actual 
line and its proposed models is clearly desirable. If a t.l.m. 
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could be constructed to exactly duplicate in magnitude and phase the 
impedance function of the actual line for all frequencies, then the 
model could be said to have "perfect accuracy". Perfect accuracy 
would probably require more lumped element sections than economically 
practical, so that some compromise is required. 

As the first, and perhaps the simplest, approach toward 
answering the questions above, this design uses as a comparison 
"device" the poles of the driving point impedance functions. The 
poles of a given length of actual line are compared to the poles of 
a proposed m-section model. The number of sections, m, in the model 
is "adjusted" until the frequency difference between corresponding 
poles of the model and actual line is sufficiently small. In each 
case, the actual line and model compared are the transposed, 
symmetrical -component equivalents of the three-phase systems. This 
simplification is sufficient for design analysis since it is not the 
detailed nature (e.g. "splitting") of the poles that is sought, but 
rather, a reasonably accurate estimate of their mean position in the 
frequency spectrum. Tlje mean positions can then be compared to 
estimate the difference between the actual line and model. The 
transposed, symmetrical -component approximation is applied to both 
the matnematical model of the actual line and that of the t.l.m., 
so that errors introduced by this approximation are self compensating. 

In analyzing the driving point impedance function of both 
the line and models, it is convenient and sufficient to consider 
two possible extremes as terminations. Both the sources and loads 
connected to the terminals of the line are non-linear devices (at 
least under some conditions of operation or fault) so that, as 
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viewed from the line, the terminals of the line may simultaneously 
appear as a high impedance to some transient components and a low 
impedance to other components. This "high-low" terminal impedance 
characteristic is provided for in the analysis to follow by consider- 
ing both open and short circuit terminations. The decision as to 
what is a "sufficiently" small difference betv/een the line and the 
model is deferred to a point later in the analysis where it can be 
made to appear a little less arbitrary. 

3.2 Driving Point Impedance of Continuous Lines 



continuous lines are well known and will not be derived for this 
discussion. Readers interested in the details of the derivation 
may consult any good text on transmission lines such as references 
7, 8, or 9. For the general transmission line circuit shown in 
Figure 3.2, the pertinent impedance relations for a length, £, are: 



The relations describing the driving point impedance of 



Z. ^ tanh yl 

s . c. 0 ' 

Z ^ = Z^ coth Y^'* 

o.c. 0 ' 



(3-1) 

(3-2) 



>Jnere 



Zc - = Impedance at Z = £ with a short 

S • C • 

circuit termination. 



Z = Impedance at Z = £ with an open 

0 • w • 

circuit termination 




(3-3) 



Y = /(R + sL)sC 



(3-4) 
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s = Complex piiase constant 

R = Positive or zero-sequence resistance per unit length 
L = Positive or zero-sequence inductance per unit length 
C = Positive or zero-sequence capacitance per unit length 
Tne shunt conductance, G, has been ignored in these relations. In 
trie short circuit case, tne poles are seen to occur when 

cosh = 0* (3-5) 

Setting y = JY-j > equation (3-5) can be expanded as follov/s: 
cosh(Y^ + JY^.)il = 0 

cosh(Y^^-)cosh(jY^A) + sinh(Yj,5.)sinh(jY^il) = 0 
cosh(YpJi)cos Y^^- j sinh(Y^Jl)sin Yj2- = 0- 
Thus, for cosh y^ = 0 it is necessary that 

cosh (Yj,Ji) cos Y-j^- = 0* (3-6) 

and 

sinh(YpJi)sin y-j^ = (3-7) 

Case No. 1 : 

Suppose Yp = 0 

Then eq. (3-7) is satisfied and eq. (3-6) 

implies y^^ = JY^^ = j(2n+l )ti/2 n=0,l,2 (3-8) 

Case No. 2; 

Suppose Yp 0 but Y^^- = i. ^ n=0,l,2 

Then eqs. (3-6) and (3-7) are never both satisfied. 

Therefore Case No. 1 gives the only solutions for y which result 
in the poles of equation (3-1). A similar development apolied to 
equation (3-2) gives for the open line: 
yJI = + jmi n = 0,1 ,2. . . . 



/ 



(3-9) 
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If equations (3-8) and (3-9) are expanded and solved for the complex 
piiase constant, s, there results: 



s.c. 






R (2n+l)^ TT^ 



ZU 



LC 



o.c. 









2U 



2 

r LC 



n=0,l ,2. . . . 
n=0,l ,2 



If tile line shown in Figure 3.2 is lossless (R = 0) then: 



^s.c. ' J^o 
^o.c. ° 




B = (jo/lT = pliase constant 



(3-10) 

(3-11) 

(3-12) 

(3-13) 

(3-14) 

(3-15) 



- . (2n + 1 )tt 
0) = + ^ — 

~ Zl AZ 

u) = + 

~ I Ac 



n=0,l ,2. . . . 
n=0,l ,2 



(3-16) 

(3-17) 



3.3 Driving Point Impedance of L-Section Models 

A schematic diagram of an L-section model for modeling the 
line of Figure 3.2 is shown in Figure 3.3. Both figures refer to 
the symmetrical component equivalent circuits of their respective 
three-phase systems. Thus, both figures represent a "transposed- 
line-approximation" to their respective non- transposed, three-phase 
counterparts. In the case that the L-section line is shorted, a 
clear distinction must be made as to which end is the source and 
wnich is the termination. In the first case, a short can be 
placed across tne last capacitor-removing it from the circuit. This 
case is referred to as "inductor input". In the second case, the 
short can be placed to connect the "last" inductor to ground. This 
case is referred to as "capacitor input". 
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3.3.1 Inductor Input L-Section ilodel (Shorted) 

The inductor input case is that shown in Figure 3.3. The 
derivation of the driving point impedance expression is rather 
lengthy and is shown in detail in Appendix B.l. The lumped element 
parameters of Figure 3.3 are defined in terms of the distributed 
parameters of Figure 3.2 as: 



R = ^ 
s m 



(3-18) 
(3-19) 
(3-20) 

Where "m" is the number of sections in the model. The result from 
Appendix B.l for an m-section; lossy model is: 



L 

s m 

P _ a 
^s ■ m • 



(R-+L )sinh mX 

7 - ^ ^ 

s.c. ~ sinh mX-sinh(m-l )X * 

The denominator of equation (B-16) is a polynomial in the real 



(B-16) 



variable cosh X so that the solutions, [cosh to the equation 

sinh mX-sinh(m-l )X = 0 yield the propagation constants in the 
propagation relation 



sR C 



(B-7) 



Solving equation (B-7) for s gives: 



s.c. ■ - (ju) i V fe' - ^)rootJ 

If the model of Figure 3.3 is lossless (R^ = 0), 

jojL^ sinh mX 

^s.c. " sinh mX-sinh(m-l )X 



m 1,2)3**** 



(B-17) 



33 




(B-18) 




(B-27) 



For an example application, suppose their are seven sections 
in a proposed model. Then m=7; there are six values of [cosh ^]pQQ^» 
and six values of s or w . The six poles of the model given by 

equations (B-26) or (B-27) are compared to the first six poles of the 
actual line given by equations (3-10) or (3-16). The difference in 
frequency between two corresponding poles is a measure of the error 
of the model for that pole number. Equations (B-16), (B-26), and 
(3-10) were implemented on the computer, as shown in Appendices C and 
D, for various numbers of sections up to about forty. A plot of the 
resulting pole error vs. pole number is shown in Figure 3.4a. The 
percent error in Figure 3.4a is defined as 



The parameters used in equations (B-26) and (3-10) for Figure 3.4a 
are the positive sequence values of a 765 kV "typical aluminum 
tower" configuration: 



% Error = 100 



I^[equation (B-26)-I^[equation (3-10)] 



(3-21) 



Ij^Lequation (3-10)J 



C = .021 pf/mi 



Z = 100 miles 



L = 1.45 mh/mi . 



R = .0251 ohms /mi 



(3-22) 

(3-23) 

(3-24) 

(3-25) 



The pole number shown in Figure 3.4a is, of course, a discreet integer 
variable with pole no. 1 corresponding to the pole of lowest 
frequency in either the actual line or model. The curve for a given 
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Plot of Pole Error vs. Pole No. 
L-Section Line (Pos. Seq.) 
Inductor Input 
Shorted Termination 
100 Mile Length 
Lossv Line and Model 
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21 Sections 
39 Sections 



Figure 3.4a 
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number of sections in the model is a series of discreet points but is 
drawn continuous to avoid confusion. Figure 3.4a shows, for example, 
that at least twelve sections are required to model the line if it is 
desired to keep all of the first five poles of the model within five 
percent error. 

If the typical parameter values of equations (3-22)-(3-25) 
are substituted directly into equations (3-10) or (B-26), it is found 

D 

that the term ^ is of order 10' whereas the terms 



(2n+l)V ^ 









[1-cosh 



7 8 

are of order 10 or 10 for even the lowest values of m or n. This 
suggests that for purposes of a pole error analysis, the resistance 
terms can be neglected when considering the positive sequence of the 
lines. This, in turn suggests that a pole error analysis using the 
lossless equations (B-17), (B-27), and (3-16) should give results 
very similar to the lossy case shown in Figure 3.4a. The results 
for the lossless case paralleling the analysis leading to Figure 
3.4a are shown in Figure 3.4b. The percent error in Figure 3.4b 
is defined as 



m 



SL J LC tl-(cosh 



(2n+l )Tf 

2 £ AX 



% Error = 100 



17n~+lT?r~ 

21 AX 



(3-26) 



As expected. Figures 3.4a and 3.4b are practically identical so that 
losses will be ignored throughout the remainder of the pole error 
analysis. Note that equation (3-26) is independent of 5-, L, and C, 
so that the pole error curves of Figure 3.4b apply to any line and 
model provided the model is "inductor input", L-section, and both 
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line and model are terminated in a short. The pole frequencies 
corresponding to the first fifteen poles of a shorted, typical 
765 kV line, for various line lengths are shown in Table 3.1. 

Equations (B-26) and (3-10) for the configuration of Figure 
3.3 were also run using zero sequence values of a typical 100 mile 
765 kV line. The pole error results for this case are shown in 
Figure 3.5 with 



i = 100 miles 


(3-27) 


R = .51 ohms/mile 


(3-28) 


C = .0143 yf/mile 


(3-29) 


L = 4.2 mh/mile . 


(3-30) 



The percent error in Figure 3.5 is as defined by equation (3-21). A 
comparison of Figures 3.4a and 3.5 shows that the pole error for zero 
sequence parameters is practically identical to the error for positive 
sequence parameters. In addition. Figure 3.5 is practically identical 
with the lossless results shown in Figure 3.4b. Thus, throughout the 
remainder of the pole error analysis, the lossy line and zero sequence 
pole error will be assumed to be the same as the lossless line error 
of any given configuration. 

3.3.2 Capacitor Input L-Section Model (Shorted) 

If the model of Figure 3.3 is turned around so that the 
capacitor is at the input, the result (lossless) is shown in Figure 
3.6. The derivation of the driving point impedance of this line is 
carried out in Appendix B.2. The analysis considers only the lossless 
case since, as discussed in section 3.3.1, the effect of resistance 
can be neglected in the pole error analysis. The results from 
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Table 3.1 

Short Circuit Pole Frequencies (Pos. Seq.) 



60 Mi 


100 Mi 


760 Hz 


456 Hz 


2,280 


1,370 


3,800 


2,280 


5,345 


3,200 


6,840 


4,100 


8,335 


5,020 


9,850 


5,920 


11,390 


6,840 


12,900 


7,750 


14,400 


8,650 


15,950 


9,580 


17,410 


10,470 


18,950 


11 ,390 


20,500 


12,300 


22,000 


13,200 



200 Mi 


300 Mi 


228 Hz 


152 Hz 


685 


456 


1 ,140 


760 


1,600 


1,054 


2,050 


1 ,370 


2,560 


1 ,675 


2,960 


1 ,975 


3,420 


2,275 


3,880 


2,580 


4,325 


2,880 


4,790 


3,190 


5,235 


3,490 


5,695 


3,790 


6,150 


4,100 


6,600 


4,400 



Table of Pole Frequencies for Shorted, Typical 



400 Mi 
114 Hz 
343 
570 
800 
1 ,025 
1 ,257 
1 ,480 
1 ,710 
1 ,940 
2,160 
2,390 
2,620 
2,840 
3,075 
3,300 



765 kV Line (Positive Sequence-Lossless) 
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Figure 3.5 
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Appendix B.2 are. 



^s.c. sinh(m+l )A-sinh mX "’“1 >2 



jojL^ sinh mX 



(B-39) 




(B-33) 



s .c. 



I\/d: n-(cosh X) 




(B-40) 



Note that equation (B-39) is not the same as equation (B-17), so that 
the values of (cosh ^)pQQ^ obtained from (B-39) for use in (B-40) will 
be different from the roots obtained from (B-17) for use in (B-27). 
Equations (B-39), (B-40), and (3-16) were solved with the computer 
as shown in Appendices C and D. A plot of the resulting pole error 
vs. pole number is shown in Figure 3.7. Having neglected resistance, 
the poles of the model are given by equation (B-40) and those of 
the actual line by equation (3-16). The percent error in Figure 3.7 
is defined by equation (3-26). Figure 3.7 shows that the error 
behavior for the capacitor input case is radically different from 
that of the inductor input. The capacitor input model consistently 
under-estimates the actual pole frequency, whereas the inductor input 
model over-estimates the frequency of the lower order poles, and, 
throughout its range, yields pole frequencies higher than the 
capacitor input model. The inductor input model gives higher pole 
frequencies due to the fact that the shorted termination removes 
the last capacitor from the model -thus reducing overall model capacity 
and increasing the model's natural frequencies. Note also, that the 
capacitor input model errors are larger at lower order poles and 
smaller at higher order poles than corresponding lengths of inductor 
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input model . 

3.3.3 L-Section Model, Open Termination 

In the case of an open termination, the poles are unaffected 
by v;hich terminals are considered to be the input or output, and only 
one analysis is required. The schematic of the open circuit case is 
just Figure 3.3 with the short removed. The driving point impedance 
for this case is derived in Appendix B.3. The results are (R^ = 0): 



Where, again, use has been made of the fact that. 



m = number of sections in the model, 
a = length of the line being modeled, 
and L and C are the electrical 
parameters per unit length of any 
given line. 



Equations (B-46), (B-47), and (3-17) were solved with the aid of a 
computer as discussed in Appendices C and D. The results are 
shown in Figure 3.8 with the pole error defined as 




sinh(m+l )A-sinh mX 



(B-46) 



(B-18) 




(B-47) 




% Error = 100 




(3-31) 






As in equation (3-26), (3-31) is independent of I, L, and C, and 
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Figure 3.8 
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depends only upon line configuration. Again, the values of the roots 
are different from what they were for either of the preceeding cases. 
The pole frequencies corresponding to the first fifteen poles of an 
open, typical 765 kV line, for various line lengths is shown in 
Table 3.2. 

3.4 Driving Point Impedance of ir-Section Models 

Referring to the line shown in Figure 3.6, it is seen that 
the L-section model can be converted to a n model by a simple 
adjustment of the first and last sections. If the capacitor at 
the mth node is replaced by one of value C^/2, and another capacitor 
of value C^/2 is added at the "0"th node, the result is a ir-section 
model. A similar type of adjustment could be made using inductors 
of value Lg/2, and the result would be a T-section model. The it 
and T-models are reciprocal networks, whereas the L-model is not. 

A pole error analysis for ir-sections follows in sections 3.4.1 and 
3.4.2. The analysis for T-sections is not included since it is 
assumed that it will add little to all the cases considered. 

3.4.1 TT-Section Model (Shorted) 

The derivation of the driving point impedance relations for 
the shorted, 7i-model are carried out in Appendix B.4 and the solution 
is carried out in Appendices C and D. The result for the shorted 
TT-model is shown in Figure 3.9. The percent error in Figure 3.9 is 
defined as in equation (3-26), but again it should be noted that 
the values of [cosh differ for each line configuration and 
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Table 3.2 



Open Circuit Pole Frequencies 



60 Hi 


100 Mi 


200 Mi 


300 Mi 


400 ^ 


0 Hz 


0 Hz 


0 Hz 


0 Hz 


0 


1,510 


905 


452 


315 


226 


3,020 


1 ,810 


905 


604 


402 


4,540 


2,720 


1 ,360 


906 


680 


6,040 


3,620 


1,810 


1 ,205 


905 


7,550 


4,540 


2,270 


1 ,510 


1,135 


9,055 


5,440 


2,770 


1,810 


1 ,385 


10,580 


6,350 


3,175 


2,120 


1 ,585 


12,090 


7,250 


3,625 


2,420 


1 ,815 


13,600 


8,175 


4,092 


2,720 


2,046 


15,050 


9,050 


4,525 


3,070 


2,260 


16,600 


9,998 


4,999 


3,325 


2,500 


18,200 


10,900 


5,450 


3,640 


2,725 


19,650 


11 ,800 


5,900 


3,940 


2,950 


21,200 


12,700 


6,350 


4,240 


3,175 



Table of Pole Frequencies for Open, Typical 
765 kV Line (Positive Sequence-Lossless) 
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Figure 3.9 
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3.4.2 TT-Section Model (Open) 

The open circuit case is derived in Appendix B.5 and solved 
in Appendices C and D. The results are as shown in Figure 3.10. For 
Figure 3.10, the percent error is as defined in equation (3-31) with 
due recognition that the values of (cosh in this case are 

different from their values in the opened, L-section line. 

3.5 Conclusions 

1. The results shown in Figures 3.4a, 3.4b, and 3.5 indicate 
that for typical EHV lines and models, the percent error between lines 
and models is largely independent of losses in either. Thus, the 
percent error for either positive or zero-sequence is the same for 

a given number of model sections and a given type of termination. 

2. The results indicated by Figures 3.4b and 3.7 show that 
due to the non-reciprocal character of the L-section model, the error 
between model and line depends greatly upon which terminals of the 
model are considered input or output. In fact, the capacitor input 
case is seen to be considerably inferior to the inductor input case 
for low order poles. 

3. Comparison of the shorted L-model with the shorted tt- 
model shows that the L-model has a rather more complicated error 
pattern than the Ti-model . The L-model first over estimates the 
true pole position for low order poles - then under estimates for 
higher order. The 7r-model consistently under estimates the position 
of the actual pole. 

4. The open circuited tt and L-model s give practically 
identical pole error plots. (Actually, each ir-model has one 
additional pole higher than the highest pole of the corresponding 
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Plot of Pole Error vs. Pole No. 
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L-model . ) 

5. Any L-section mode! can be converted to either a tt or T 
model by a simple change in the first and last sections to obtain 
any of the pole-error patterns of proceeding sections. 

6. A suggested method of using the pole error plots is 
shown in Figure 3.11. Figure 3.11 shows, as an example, the pole 
error "curves" of a twelve section L-model as taken from Figures 3.4b 
and 3.8. A five percent error bound is suggested by the dotted lines 
at the top and bottom of the figure. The region between the solid 
curves is the "error band" of the twelve section L-model. If the 
terminals of the line appear as a predominantly high impedance to 
frequencies corresponding to the first five poles, the line error 

is given (approximately) by the "open" curve. A corresponding 
observation applies for the low impedance case. The dashed curve 
mid-way between the solid curves is meant to indicate the "nominal 
error" when the terminals appear as neither open nor shorted. If 
twelve sections were used to model a sixty mile line. Figure 3.11 
shows that the error band lies within the five percent error 
bounds up to about the fifth pole. As seen from Tables 3.1 and 3.2 
the fifth pole corresponds to about 6000 Hz for a sixty mile line. 

One could say then - "the 12-section L-model is good to within five 
percent up to 6000 Hz" (for a 765 kV line). An error plot 
similar to that of Figure 3.11 could be made for any model type 
composed of any number of sections. In particular, if it is 
desired to model a 400 mile line to "within five percent" an 
inspection of Figures 3.4b, and 3.8 shows that a model of about 
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Example of an "Error Band-Error Bound" 
Plot for a Twelve Section L-Model 



Figure 3.11 
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forty sections will ensure that the line is "good to within five 
percent" up to the fifteenth pole. Tables 3.1 and 3.2 show that 

the fifteenth pole corresponds to a frequency of about 3100 Hz 

(for a typical 765 kV line). Thus, a minimum of forty ten-mile 

sections would be required to model a 400 mile line to within the 

five percent criteria. 

Tables 3.1 and 3.2 give the pole frequencies, specifically, 
for a typical 765 kV line. Actually, as will be shown in Chapter 4, 
corresponding tables for most other EHV lines would be very similar 
to Tables 3.1 and 3.2 (within about ten percent). With Tables 3.1 
and 3.2 as representative of EHV lines, it is now possible to make 
some quantitative decisions concerning the number of sections to be 
built for the t.l .m. 

7. As stated earlier in this report - as a minimum - the t.l.m. 
is to model EHV lines from sixty to 400 miles in length. The five 
percent "specification" referred to earlier is actually a result of 
observations such as those in Conclusion 6. It was noted earlier 
that forty sections is the minimum number which will just ensure 
that a 400 mile line-model is within five percent to 3060 Hz. This 
in turn, implies at least forty ten-mile sections for the model. In 
addition, it is proposed to model all line lengths between sixty 
and 400 miles by changing the number of sections in the model. Thus, 
if a line to be modeled is 396 miles long, it would be desirable to 
set up a model "396 miles long". This adjustment could be 
accomplished with the ten-mile sections by using "taps" to adjust 
length. However, this t.l.m. design must model, horizontal or 
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vertical, non- transposed lines - this too requires taps. To avoid 

having to tap the sections for length, it is proposed to make the 

sections "fine" enough (i.e. short enough) so as to make tapping 

unnecessary for length adjustments. Line lengths to 400 miles can 

be reached to within one percent by making the largest sections 

have only eight-mile length. Thus, the 396 mile line, e.g., could 

be reached (to one percent) by either forty nine or fifty eight- 

mile sections. But eight-mile sections are not "fine" enough to 

model short lines, etc., etc. Finally then, with considerations 

such as those above, it was decided that the t.l.m. design should 

consist of thirty four-mile sections and thirty-five eight-mile 

sections. This selection is seen to be quite conservative (by 

the chosen criteria) and provides a model somewhat "finer" than 

1-4 

some other t.l .m. 's. 

8. The pole error analysis is seen to provide a means of 
quantitatively comparing proposed models, and is useful in indicat- 
ing the relative merit of, say, doubling the number of sections - 
or halving them, etc. 

3.6 Applications and Recommendations 

The solution technique employed to solve the impedance 
functions of this chapter can have application in other fields of 
study. The denominator of each of the impedance functions of this 
chapter corresponds to the "frequency polynomial" of mechanical 
vibrations problems. If the frequency polynomial is set equal to 
zero, it becomes the "frequency equation" for determining the 
natural frequencies of mechanical systems. Even in systems with 
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only two degrees of freedom the frequency equation is of fourth degree 
(quadratic in w ) so that, in continuous systems or systems having 
more than two degrees of freedom, complicated solution techniques 
must be employed or approximations made. The technique suggested 
in this chapter and employed in Appendices C and D is a means of 
having the computer generate the frequency equation and then solve 
it for the natural frequencies of the system. A detailed investiga- 
tion into the exact class of problems for which this solution technique 
may have validity has not been conducted by the author, but it is 
suggested that such an investigation could be worthwhile. 
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CHAPTER 4 

PARAMETERS OF EHV LINES AND MODELS 
4.1 Parameters of Representative EHV Lines 

In order to determine the parameters or component values for 
a t.l.m., it is first necessary to know the parameters of typical 
lines the t.l.m. would be expected to model. The t.l.m. should be 
capable of modeling not only lines already in existence but, in 
addition, it should be capable of modeling proposed and even hypothetical 
lines. For information pertaining to the parameters of existing lines 
use was made of Report No. 710, Electrical Characteristics of Trans- 
mission Lines, of the American Electric Power Service Corporation 
(AEP)^^. Appendix IV of Report No. 710 describes in detail thirty- 
four different EHV transmission line configurations in use by the 
company. From the tiiirty-four configurations a sample of six was 
selected as representative of types of configurations. For example, 

AEP type "H" was taken as one representative of a 345 kV, single 
conductor per phase, single ground wire, vertical configuration. 

Type "H-Special" was selected since it is a common line and represents 
a form of triangular configuration. A hypothetical type "39F" was 
"constructed" by using two phase-conductors per bundle spaced 
eighteen inches apart - this was taken as representative of a two- 
conductor, "vertical", 345 kV line. A type "41A" line was selected 
to represent a 345 kV, two-conductor, horizontal configuration. The 
AEP "Typical Aluminum Tower" configuration was taken as representative 
of most 765 kV lines (horizontal). A hypothetical 765 kV "H- 
Special" was constructed by scaling the 345 kV "H-Special" conductor 
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distances in the ratios that exist between the type "41 A" and the 
"Typical Aluminum Tower". Finally, two hypothetical 1500 kV con- 
figurations were taken from reference 11 - a horizontal configuration 
(Geometry No. 2 and Bundle Drawing No. 4), and a triangular line 

(Geometry No. 3 and Bundle Drawing No. 4). 

Typical positive and zero-sequence parameters obtained for 

these representative lines are shown in Tables 4.1 and 4.2. The 

values shown in Tables 4.1 and 4.2 (except for X^/X|^ and Q) were 

obtained with the aid of a computer using the "Electric Power Systems 

1 2 

Engineering Laboratory (EPSEL) Line Parameters Program" This 
program requires as input the geometrical configuration of the 
conductors of a given line (along with conductor and earth 
properties) and computes the inductance, capacitance, and resistance 
matrices of the actual line as a function of frequency. In addition, 
as well as the parameters of Tables 4.1 and 4.2, the program computes 
and prints the values of the model parameters of Figure 2.2. The 
model parameters are printed on the basis of a one-mile long section 
and assumes a scale factor of one betv/een impedances of the line and 
model. Model parameter values at sixty Hz, on the basis described 
above, for the eight representati ve line models are shown in Table 
4.3 (earth resistivity is taken to be 100 ohms-meters) . Referring 
to Table 4.2 note that, except for the 1500 kV lines, the positive- 
sequence (or average) Q at 60 Hz is about twenty or lower. The higfi 
Q values for the 1500 kV lines are questionable and would be 
difficult to duplicate in a practical model. Consequently, no 
attempt will be made to even approach the Q values suggested by the 
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Pos. Seq. 
1/Xo 



Pos. Seq. 
X, 



Pos. Seq. 
R 



Type 

39F 



7*11 0*594-7?- 0.063 - 7 ?- 236,000 (mi. 

meg si mile mile » \ / 



345 kV 
H-Spec. 



5.25 



0.810 



0.087 



263,000 



Type 

H 



5.65 



0.751 



0.087 



236,000 



765 kV 
H-Spec. 



7.49 



0.561 



0.027 



239,000 



Type 

41A 



6.82 



0.627 



0.054 



234,000 



1500 kV 
Hori zontal 



8.08 



0.522 



0.005 



237,000 



1500 kV 
Triangular 



8.41 



0.494 



0.005 



240,000 



765 kV 
Horizontal 



7.78 



0.548 



0.028 



236,000 



Table of Positive Sequence Parameters for 
Some Typical EHV Line Configurations (at 60 Hz). 



Table 4.1 



I 




jff* 
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Type 

39F 



345 kV 
H-Spec. 



Type 

H 



765 kV 
H-Spec. 



Type 

41A 



1500 kV 
Horizontal 



1500 kV 
Triangular 



765 kV 
Horizontal 



Zero Seq. 



’•59 It 

mi 



1.85 



1.97 
1.66 
1.79 
1.10 
1.23 
1 .64 



Zero Seq. 
R 



■329 ^ 
mi 



0.395 



0.395 

0.418 

0.598 

0.131 

0.133 

0.522 



Zero. Seq. 

1/Xp 



3.92 



mi 



meg n 



3.63 



3.30 



4.65 



4.80 



5.83 



4.81 



5.31 



Pos. Seq. 

Q 

9 

9 

9 

21 

12 

104 

99 

20 



Table of Positive Sequence Q and Zero Sequence 
Parameters of some Typical EHV Line 
Configurations (at 60 Hz). 



Table 4.2 
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rameter 


Type 

39F 


Ml 


1.83 


42 


1.76 


43 


1.59 


*-00 


0.736 


‘-12 


0.297 


1' 

‘-23 


0.147 


f ' 
‘'12 


3.6 


r ' 

‘'23 


3.5 


c 

‘'13 


1.5 


c 

‘'10 


11.2 


c 

‘'20 


9.2 




10.8 


"il 


56.2 


42 


60.4 


42 


71.4 


4o 


89.0 



345 kV 


Type 


H-Spec. 


H 


2.21 


2.21 


2.16 


2.16 


2.05 


2.05 


0.928 


0.928 


0.022 


0.276 


-0.051 


0.168 


1.3 


2.5 


1.5 


2.6 


1.5 


1.2 


10.0 


9.5 


9.6 


8.0 


9.7 


8.9 


80.0 


80.0 


84.6 


84.6 


94.9 


94.9 


103. 


103. 



Table of Unsealed Parameters for Models 
of Representative EHV Lines. (L‘ in mh/mi . 
C in nf/mi . , and R' in mfi/mi . ) 



765 kV 
H-Spec. 

1.52 
1 .49 
1.42 
0.975 
0.004 
-0.038 
2.2 
2.7 
2.7 

13.2 

12.2 

11.5 

16.5 
24.0 
40.4 

130. 



Table 4.3 
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Parameter 


Type 

41A 


1500 kV 
Hori zontal 


1500 kV 
Triangular 


765 kV 
Horizontal 


Ml 


1.81 


1.52 


1.30 


1.60 


M2 


1.78 


1.40 


1 .46 


1 .59 


1 ' 
^33 


1.81 


1.52 


1.30 


1.60 


'•00 


0.886 


0.410 


0.604 


0.820 


M2 


0.212 


0.147 


0.066 


0.214 


1 ' 
^23 


0.212 


0.147 


0.066 


0.214 


C 

^12 


2.4 


2.7 


3.1 


2.95 


r ' 

^23 


2.4 


2.7 


3.1 


2.95 


C ' 
^13 


0.7 


0.7 


3.3 


0.76 


Mo 


13.2 


15.8 


12.5 


14.7 


*'20 


11.9 


14.8 


13.3 


13.1 


*'30 


13.2 


15.8 


12.5 


14.7 


Ml 


51.5 


6.2 


4.7 


26.0 


R'2 


58.0 


4.0 


6.7 


30.6 


Ms 


51 .5 


6.2 


4.7 


26.0 


0 
- 0 
a: 


182. 


41.5 


42.3 


165. 




Table of Unsealed 


Parameters 


for Models 





of Representative EHV Lines. (L‘ in mh/mi . , 
C in nf/mi . , and R in mf2/mi.) 

Table 4.3 (Cont. ) 
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1500 kV line values. With a Q of twenty, the 765 kV horizontal line 
will be taken as the minimum standard (in Q) to be achieved by this 
t.l.m. This means (referring to Table 4.2) that each of the sixty- 
five sections of the t.l.m. must have a positive sequence Q of at 
least twenty at sixty Hz. The positive and zero-sequence Q of 
transmission lines is not constant with frequency and, in fact, 
increases with frequency. Graphs of the positive and zero-sequence 
Q of representative lines are shown in Figures 4.1 and 4.2. The 
positive-sequence Q is seen to increase steadily in the range to 
3060 Hz while the zero-sequence Q becomes nominally constant in the 
range from 1500-3060 Hz. Each model section must provide a Q 
characteristic equivalent to those of Figures 4.1 and 4.2 in order 
to accurately model the actual line. If the model section is 
designed to meet the positive-sequence Q requirement of the typical 
765 kV line, than it can be made to model lines of lower Q by 
simply inserting resistance in the phase conductors. 

The positive-sequence inductance of EHV transmission lines 
is practically constant over the frequency range of interest, but 
the zero-sequence inductance is not. A graph of the variation of 
zero sequence inductance with frequency is shown in Figure 4.3 for 
the typical 765 kV line. The zero-sequence inductance is seen to 
drop markedly with increasing frequency-parti cularly in the low 
frequency region. Thus, in the model also, the zero-sequence 
inductance must correctly decrease with increasing frequency and 
produce (as nearly as possible) the scaled equivalent of Figure 4.3. 



Graph of Positive Sequence Q vs. Frequency for Representative EHV Lines 
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* 

a. 

cr 



Figure 4.1 



I 






Graph of Zero-Sequence Q vs. Frequency for Representative EfIV Lines 
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Figure 4.2 




1 



liraph of Zero-Sequence Inductance vs. Frequency for Typical 765 kV Line 
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Figure 4.3 
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4.2J Line/Model Scale Factors and General Model Parameters 

Due to the fact that the t.I.m. must be compatible with existing 
TSS transformer and compensating reactor models, the base voltage of the 
t.I.m. must be 138 V line-to-line, and the ratio of the base impedance 
of the model to that of a typical 765 kV line must be 0.355. This means 
that the base current and power of the t.I.m. are 38.6 ma and 9.2 watts 
respectively. The base impedance of a typical 765 kV line on a 100 MVA 
base is 5850 ohms. The base impedance of the model is 2078 ohms. With 
the above impedance ratio constraint between the model and 765 kV line. 



and 



*•8 model 
**B actual 

model 
*^B actual 



.355, 

.355, 



B model _ 1_ 

^B actual 



= 2.82 . 



To model the typical 765 kV line, all the 765 kV parameters of Table 
4.3 must be appropriately scaled by the ratios above. For example, 
L 22 of Table 4.3 becomes .564 mh/mi , or 2.26 mh/4 mi, or 4.52 mh/8 mi 
Thus, the 4-mile section would be constructed with a basic 2.26 mh 
inductor. An identical procedure applied to the 345 kV type "41A" 
line gives (on a 100 MVA base): 



^^B^345 ^ H90f2 

2078 . , ,, 

- 1190- 

L 22 =3.1 mh/mi, or 12.4 mh/4 mi, etc. 

To model the 765 kV line the inductor in the four mile section must 
be 2.26 mh and to model the 345 kV line it must be 12.4 mh. This 
implies a very wide range of taps on the model inductors. A wide 



i 



A 



4 
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range of taps on an inductor is not only difficult and expensive to 
construct, but has the further disadvantage that the Q of the inductor 
gets "tapped" as the inductance. Thus, tapping the inductance down to 
one-sixth of the total inductance cuts the Q down to about one-sixth 
of what it was for the total inductor. 

A procedure for solving this problem is to choose, for lines 
other than the typical 765 kV, a power base other than the commonly 
used 100 HVA. Again suppose that the 345 kV type "41 A" line is to be 
modeled. The first step is to establish an impedance ratio between 
the actual line and model such that 1^2 ^or the 345 kV model is the 
same as that already provided for the typical (horizontal) 765 kV 
line. In other v/ords, using Table 4.3, seek a scale factor F, 



c _ model 

■ line 



(4-1) 



such that 



''^*-^^345 ' ^*-22^ 765 hor. 



(4-2) 



So that 




(4-3) 



and 



But 



so that 




(4-4) 



(4-5) 



Thus, when modeling a 345 kV type "41A" line, it is necessary to use 
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a power base in the line of 18.1 MVA instead of the usual 100 MVA. 

The resulting parameters required of a four-mile model section are 
shown in Tables 4.4a and 4.4b respectively. Note in Table 4.4a that 
the inductors corresponding to phases "a" and "c" are equal - reflect- 
ing the symmetry of the horizontal lines. This is not the case for 
the "vertical" lines of Table 4.4b. With the inductances of the 765 
kV horizontal model as a base, note that all the L^'s of Table 4.4b 
are larger than 2.28 and all the L^'s are smaller. In fact, all the 
L,'s of Table 4.4b are larger than 2.28 by the same amount (A) that 
all the L^'s of Table 4.4b are smaller than 2.28. The scale factor 
to achieve this "balance" for vertical configurations is not given 
by equation (4-3) but is determined as follows: 

Using Table 4.3, seek a scale factor F such that 



The use of equation (4-6) ensures that the unbalance inherent in 
vertical lines is shared equally by both phases "a" and "b" in the 
model. A "balanced sharing" is desirable since a factor determined 
as in equation (4-3) gives the result that the phase "c" model 
inductor would be required to span a larger range than the phase "a" 
inductor - thus requiring more taps. Use of equation (4-6) ensures 
that the phase "a" and "c" inductors will have the same number of 
taps. Referring again to Table 4.4a, note that for the two 




- A. 



Solving for F [using the fact that (L^,) 



11^765 hor." ^^33^765 hor.^’ 



F = ^^^11^765 hor. 
^^11 '43^ vert. 



(4-6) 



f 
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Parameter 


Type 

41A 


1500 kV 1500 kV 

Horizontal Triangular 


765 kV 
Horizontal 


'•a 


2.29 


2.34 


2.16 


2.28 


Lb 


2.26 


2.16 


2.41 


2.26 


Lc 


2.29 


2.34 


2.16 


2.28 


Ln 


1.12 


0.630 


1 .0 


1.16 


’■ab 


0.27 


0.23 


0.11 


0.30 


Lbc 


0.27 


0.23 


0.11 


0. 30 


Lab 


30.4 


28.1 


30 


33.2 


Sc 


30.4 


28.1 


30 


33.2 


Sc 


8.8 


7.0 


32 


8.6 


^an 


166 


164 


121 


166 


‘'bn 


149 


154 


129 


148 




166 


164 


121 


166 


^'a 


65 


10 


7.85 


37 


«b 


73 


6.5 


10.8 


43 




65 


10 


7.85 


37 


R 

n 


230 


67 


70 


235 




Table of Scaled Parameters 
4-i1ile Model Section. (L 
C in nf/4 mi . , R in 


Implied for a 
in mh/4 mi . , 
mSV4 mi . ) 





Table 4.4a 
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Parameter 


Type 


345 kV 


Type 


765 kV 




39 


H-Spec. 


H 


H-Spec. 


'a 


2.44 


2.36 


2.36 


2.34 


l-b 


2.34 


2.31 


2.31 


2.29 


Lc 


2.12 


2.20 


2.20 


2.2 


'■n 


0.98 


0.99 


0.99 


1.50 


*'ab 


0.40 


0.24 


0.30 


0.006 


Sc 


0.20 


-0.05 


0.18 


-0.06 


‘'ab 


43.2 


19.4 


37.3 


22.8 


Sc 


42.0 


22.4 


38.8 


28.1 


o 

o 


18.0 


22.4 


17.9 


28.1 


Sn 


135 


149 


142 


137 


^bm 


111 


143 


119 


127 


'cn 


129 


145 


131 


119 




75 


86 


86 


20 


«b 


80 


91 


91 


37 




95 


102 


102 


62 


f'n 


117 


no 


no 


201 


Table 


of Scaled 


Parameters Implied 


for a 4-Mile 





Model Section. (L in mh/4 mi., C in nf/4 mi., R in mQ/4 mi.) 



Table 4.4b 
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1500 kV lines is not 2.26 mh as would be found had equation (4-3) been 
used to determine the scale factor. In fact, still another technique 
has been employed to determine the scale factor of these two lines. 

This third technique will be discussed later in this chapter. 

Some remarks are in order at this point concerning the model 
zero-sequence inductance, ^ , and the inductor of Table 4.4a. 

To obtain the zero-sequence characteristic required of a 4-mile model 
section for modeling a 765 kV line, for example, the graph of Figure 
4.3 is simply scaled by the factor (0.355)(4.0) = 1.42; the result 
is shown in Figure 4.4. The model zero-sequence inductance is required 
to fall from 6.2 mh at 60 Hz to about 4.0 mh at 1000 Hz. The point 
to note is that of Tables 4.4a-b is not the zero-sequence inductance 
of the model per se, but is the component which must be designed 
to produce the characteristic shown in Figure 4.4 if a transposed 
765 kV horizontal line is to be modeled. The detailed design of a 
suitable for this t.l.m. is being carried out by Hr. William Feero 
as part of his thesis to be completed in September, 1971. Thus, the 
design of is being carried out in parallel with (but slightly 
behind) the design of the remainder of the t.l.m. 

In preparation for the design of the model inductors, 
excluding consideration of the 1500 kV lines of Tables 4.4a and 4.4b, 
it was decided that the inductance values shown in Table 4.5 should 
be provided in the 4-mile sections in order to model most EHV 
lines. 






> 




Modeling a Typical 765 kV Line 
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Figure 4.4 
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Table 4.5 

Self Inductance Designed Into the 4-Mile Sections 



!:a 


hb 




2.24 mh 


2.16 mh 


2.08 mh 


2.29 


2.21 


2.13 


2.34 


2.26 


2.18 


2.39 


2.31 


2.23 


2.44 


2.36 


2.28 


2.49 


2.41 


2.33 



Early in the design of the t.l.m. the 1500 kV lines were 
generally ignored as beyond the modeling capabilities attainable by 
this design. With respect to the quality factor, the preceeding 
statement is, indeed, true. However, in all other respects it is 
possible to model these lines also-with a proper choice of scale 
factor. Note in Table 4.5, under L|^, the value 2.26 mh originally 
provided for modeling the typical 765 kV line. Equation (4-3) shows 
the L 22 (unsealed Lj^) of a 345 kV line being scaled to the 2.26 mh 
inductor previously provided. The 2.26 mh value is not unique, and 
any inductance parameter of a line to be modeled could be scaled 
to any one of the inductances provided in Table 4.5. In fact, to 
properly fit the two 1500 kV lines, the L|^ of the horizontal 1500 
kV line was scaled to 2.16 mh and the triangular to 2.41 mh. 

(The extremes of the existing L|^ were used in order to keep the 
resulting values of (Lj.)-| 5 Qg and (Lg)i 5 gg within "reach" of the 
ranges provided in Table 4.5.) The above procedure resulted in 
an impedance scale ratio of 0.414 for the triangular line and 0.384 
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for the horizontal. The resulting power bases to be applied to the 
actual lines are 448 MVA for the triangular, and 417 MVA for the 
horizontal. Thus, the three scaling examples provided in this 
section illustrate the considerable flexibility of the "Floating 
Power Base" technique. The values of self inductance to be designed 
into the 8-mile sections are obtained by simply doubling the values 
in Table 4.5. 

4.2.2 Selection of Hodel Mutual Inductance Parameters 

Tables 4.4a and 4.4b indicate that typical values of 
range from 0.006 mh/4 mi to 0.4 mh/4 mi, and those for from 
-0.06 mh/4 mi to 0.27 mh/4 mi. Providing this wide range of values 
is neither practical nor necessary for the following reasons: 

1. To cover the range implied above would require too 
many taps on the mutual windings and too many values 
brought to the front panel of the section. As few 

as six "tap values" per component in Table 4.4 means 
that 120 lines must be brought to the front panel. 

(It is proposed to limit the number of taps on any 
component to six.) 

2. In those cases where the mutual inductance value 
is 1/20 or less of the self inductance upon which it 
is wound, it would be nearly impossible to measure 
the effect of the mutual inductance. Therefore, 
those values of mutual inductance amounting to 1/20 
or less of the self inductance in the same line will 
simply be ignored. Note that this rule leaves the 
negative mutuals ignored. 
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3. Discounting the extremes mentioned above, most 
mutual values appear to be in the range of from 
0.17 mh to about 0.27 mh - the single exception 
being the far removed 0.40 mh of the hypothetical 
type "39F" line. Mutual taps actually provided in 
the model section should reflect this "crowding" in 
the mid-band by the provision of relatively more 
taps in this band. 

4. Horizontal lines are common, and it is seen that 
for horizontal lines the mutual inductances are equal. 

The selection of taps to be provided in the model 
section should preserve this balance if possible. 

Based upon the arguments above, the values of mutual inductance 
to be designed into the 4-mile model sections are shown in Table 
4.6 - those for the 8-mile sections are shown in Table 5.5. 

Table 4.6 





L^^(rah) 


cn 


.15 


CO 


CO 

• 


• 

ro 

o 


• 

ro 

o 


.23 


.23 


• 

ro 

cn 


.25 


.33 


.32 



Table of Mutual Inductance Values 
Designed into the 4-Mile Sections 
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4.2.3 Selection of Model Resistance Parameters 

The resistance ranges implied for a 4-mile model section are 
shown in Tables 4.4a and 4.4b in which the figures represent the 
total resistance required of a section from input terminal to output 
terminal. Thus, the resistance indicated includes the resistance of 
the self and mutual inductance windings discussed above. The amount 
provided for operator adjustment is the difference between winding 
and connection resistance and the values indicated in the tables. 

The resistance of the self and mutual windings must obviously be less 
than the resistance values in Tables 4.4 if the 4-mile model section 
is to provide a correct modeling of line Q. The exact value of the 
resistances to be designed into the 4-mile section cannot be deter- 
mined until the resistance of the self and mutual windings is known. 
The setting of the operator adjusted resistance parameters to be 
provided is deferred to Chapter 7. 

4.2.4 Selection of Model Capacitance Parameters 

The capacitance parameters implied for a 4-mile section are 
also indicated in Tables 4.4a and 4.4b. With a wide range of 
capacitor sizes readily available from manufacturers, a bank of 
six capacitors will be provided for each capacitance parameter for 
a total of 36 capacitors in each section. The capacitors of each 
bank are to be sized so as to span the ranges indicated with "room" 
on either side. With six capacitors per capacitance parameter, 
about 41 values can be presented to cover the range of the parameter. 
Tables of the actual capacitor banks purchased for each parameter 
are presented in Chapter 6. 
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4.3 Summary Remarks 

1. The parameters of the six non-1500 kV lines indicated 
in Tables 4.1, 4.2, and 4.3 are taken to "span" the range 
of typical EHV lines (including 500 kV lines). 

2. No attempt is made to model positive sequence line 
Q's greater than 20. 

3. Six model values spanning the range of EHV line 
configuration "types" are provided for each parameter. 

4. Only three scaling techniques were illustrated, but 
any line to be modeled could be scaled to any one of the 
existing model values in an effort to make diverse lines 
fit the model. 

5. A scale factor, F, obtained for inductance and 
resistance scaling is used as 1/F for capacitance 
scaling. 

6. Values of mutual inductance less than five percent 
of the self inductance in the same line are ignored. 

7. Parameters indicated in Tables 4.4a and 4.4b for a 
4-mile model section are simply doubled to obtain the 
values for an 8-mile section. 

8. The details of the return path inductor design, L^, 
will appear in a subsequent report (by a different author) 
due in September, 1971. 







$ 
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CHAPTER 5 
INDUCTOR DESIGN 

b,1 Model Inductor Requirements 

The model inductance values required of a 4-mile section are 
shown in Tables 4.5 and 4.6. The inductance values for an eight-mile 
section are simply double those of Tables 4.5 and 4.6. In addition 
to simply providing the inductance, the inductors in the model sections 
must have considerably greater than the Q characteristic for the 
"best" line to be modeled. The Q characteristic of the 765 kV 
Horizontal line shown in Figure 4.1 is taken as the minimum standard 
to be met by the model sections. The Q characteristics of the 
individual model inductances must exceed the requirements of the over- 
all section. While the line and section require a Q of twenty at 
60 Hz, the inductors should have a Q of perhaps 40 at 60Hz - and a 
similar excess over the entire frequency range. The characteristic 
shown in Figure 4.1 is the "average" Q of the three phases so that 
it is not the exact characteristic to be tracked by any of the 
three self inductors of the model even if this tracking were possible. 
Tiie solution of the "Q Tracking Problem" in this design is to design or 
select an inductor whose Q characteristic has the general shape shown 
in Figure 4.1 but whose Q value is approximately fifteen to twenty- 
five points higher along the range. Any excess Q of the model can 
be handled by simply inserting an additional small resistance in 
series witn the model inductors. 

In addition to the Q characteristic per se, the model 
inductors must produce the required Q at the very low values of base 
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current selected for the model, and must have an inductance that remains 
essentially constant up to at least 25 p.u. fault current. The 
inductance of all the piiase and mutual inductors must also remain 
constant with frequency - only the return path inductor, L^, is to 
vary with frequency to provide the zero sequence inductance characteristic. 
5.2 Inductor Core Selection and Core Properties 

Tne inductor characteristics described in section 5.1 are 
determined, to a very large extent, by the core material of the inductor. 
Tne tneory leading to tne final relations describing inductor behavior 
is rather involved and will not be included in this report - only the 
more significant relations will oe presented to indicate the overall 
procedure. Readers interested in a more complete background leading to 
tne relations presented herein should consult references 13-19. 

A number of core materials and configurations v/ere investigated 
to meet the specifications described above, but most were readily found 
to be unsuitable and rejected early in the design. A potted ferrite 
core, for example, meeting all other requirements, was rejected as too 
expensive (Ferroxcube series 6656). Two alternatives were considered 
in some detail - oriented silicon steel "C" cores and molybdenum- 
nickel-iron powdered toroids. The oriented silicon was rejected when 
the model base current was finally set at the relatively low value of 
38.6 rna. At this low current the oriented silicon could not produce 
tne desired Q characteristic at a sufficiently low cost or weight. 
Basically, the Q of oriented Fe-Si Inductors increases with increasing 
current (up to saturation), v;hereas with molybdenum-nickel, Q 
decreases with increasing current. Thus, the Mo-Ni cores nave their 



highest Q's at low currents and the Fe-Si cores have their lowest Q's 



79 



at low currents. The basic relation governing the Q of an inductor can 
be written as (for the frequencies of interest here) 



Q = 



toLl‘ 



R. + R 
dc ac 



Where 



(5-1) 



L = Inductance in henries 
oj = 2 it X frequency 

Rdc = dc resistance of the winding (J2) 

R^„ = equivalent core loss resistance (fi) 

gC 

2 

Witi) these definitions it is seen that R,^I corresponds to both the 

dC — 

hysteresis and eddy current losses (See reference 17 p.p. 58ff , ) and 
thus represents the total core loss. R,_, the core loss coefficient, 

aC 

1 8 

is itself a function of current. For powdered cores the manufacturer 
provides data on R explicitly, whereas for Fe-Si the entire term 

dC 

2 

R^^I is provided. If the data provided by the manufacturer is 
converted to the terms of equation (5-1), it is found that for Fe-Si 
the core loss, W, is given by (approximately) 

W = k^(w)I^•^ (5-2) 

whereas for the powdered Mo-Ui the loss is given by (approximately) 

W = R^^I^ = k2((o)I^*^ . (5-3) 

In equations (5-2) and (5-3), k^(w) and k 2 (co) are core loss constants 
and are a function of frequency. If equations (5-2) and (5-3) are 
suustituted in equation (5-1) there results, respectively. 



Q 



Fe-Si 



0 ) 11 ^ 



(5-4) 



and 
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(5-5) 






Equations (5-4) and (5-5) apply only in the region below saturation 
and suggest the reason wily the Q of the Fe-Si core increases with 
current, whereas in the Mo-Ni core Q decreases. The Q calculations 
were conducted for a number of oriented Fe-Si cores (Arnold-AH-393, 
AH-34, AL-59, AH-95, AH-86, AH-31), maximized with respect to air gap, 
with the result that cores up to eight pounds each, costing up to 
$39.00 each, could not produce the required Q characteristic at the 
base current levels required. 



smallest (in terms of those offered by the Arnold Engineering Co. or 
I'lagnetics Corp.) and least expensive unit which could be made to 
produce the desired characteristic with a Q of 40 at 60 Hz was Arnold 
part no. A-218259-2. Tiie Q characteristic of this small (2.25 inches 
O.D.) toroid when full wound is shown in Figure 5.1. The curve above 
100 Hz is the manufacturers specification which was verified for a 
typical full wound core using a GR 1650A inductance bridge down to 
300 Hz. The curve below 100 Hz is the extrapolation of the author. 
The Q at 59 Hz was checked using the circuit shown in Figure 5.2. 

was determined by measuring the dissipation factor, D, of the 
870 yf capacitor bank and using the relation 



Several Mo-Hi cores could meet the characteristic, but the 




(5-6) 



D was measured at 60 Hz using a GR 1611 -A capacitor bridge and found 
to be .0045. Five copper lug connections are included in the tank 
circuit along with (hirteen feet of No. 10 AWG solid copper magnet 



Quality Factor vs. Frequency for T.L.M. Inductors and 4-Mile Sections 
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Figure 5.1 



CIrcuU for Measuring T.L.M. Inductor Q at 59 Hz 
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wire for a total measured resistance of 15 mfi using a GR 1650-A 
resistance bridge. The voltage waveform across the inductors was 
observed with a Tektronix 549 storage oscilloscope. The 8.4 mh 
inductor shown in the figure is four 2.1 mh inductors wound on the 
t.l.m. cores. The tank circuit was "pumped up" using the oscillator 
at 59 Hz after which the switch was opened and the decay waveform 
observed on the storage oscilloscope. Twelve separate measurements 
were made covering the range of initial voltages from 0.02 v peak 
to 0.4 V peak - corresponding to peak initial currents in the range 
from 7 ma to 129 ma. In each case the time was measured for the 
voltage across the tank to fall tol^ of its initial value. The 
times measured ranged from a low of 141 m sec at 0.4 v peak to 
157 m sec below 0.1 v volt peak. The total resistance in the tank 
was determined from 



R 

•^T time 



16.8 rnh 
time 



Hhere 



R_ = R^ + R^ + R , + R 

T c H dc ac 



Finally, Q of the inductors was determined from 



(5-7) 

(5-8) 



Q = 



cjL 

^dc " "ac 



(5-9) 



The twelve Q values at 59 Hz ranged from a low of 34.6 at .4 v peak 
(initial I = 129 ma), to a high of 40 for initial voltage values 
below 0.1 V peak (32 ma). A similar "pumping" technique was used 
to measure the inter-winding capacitance of a single t.l.m. inductor 
of 2.1 mh. A series circuit was formed of an H.P. 4204A variable 
frequency oscillator, a 4.7 meg Q resistor, and a single 2.1 mh 
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t.l.m. inductor. The voltage across the coil was observed while the 
frequency was increased. The voltage peaked sharply at 300.0 KHz 
implying an inter-winding capacitance of 134 p.f. 

The maximum inductor Q shown in Figure 5.1 occurs when the 

losses in the core of the inductor equal the losses in the resistance 

of the windings. At low frequencies the winding losses predominate; 

at maximum Q the two losses are equal, and for frequencies above 

Q the core losses predominate. The initial slope of the inductor 
mdx • 

Q curve is determined by winding losses which remain essentially 
constant in this frequency range. In other words - at low frequencies, 
the Q is determined by winding resistance alone. In fact, a t.l.m. 
inductor of 2.1 mh at 60 Hz has toL = 800 mf^ and a winding resistance 
of 20-mQ indicating a Q of 40. To bring the Q of the t.l.m. section 
down to 20, for example, 20 mil of resistance is to be inserted into 
the inductor winding - sufficient to bring the overall Q of the 
section at 60 Hz down to about 20. Other, larger values of resistance 
would have to be added to produce Q's less than 20 for modeling lower 
Q lines. The second curve in Figure 5.1 represents a typical Q 
characteristic of a section . Note that the section characteristic 
is lower, and the maximum is shifted to the right - reflecting the 
fact that core losses must reach a higher value before they can 
equal the winding losses thus producing the Q_ ^ condition. Note 
also that having shown a case in which the winding resistance was 
doubled, the initial slope of the section Q curve is about one-half 
of that of the inductor alone. The third curve in Figure 5.1 is that 
of the typical 765 kV line transferred from Figure 4.1; this curve 










'C:> 
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represents the highest Q which must be achieved by a section. It is 
seen that setting the Q of the section at 60 Hz to twenty results 
in the section Q falling below that of the line above 1700 Hz. Thus, 
if Q is a critical factor above 1700 Hz in a particular test on a 
typical 765 kV line, it will be necessary to accept higher Q at sixty 
Hz by setting less than 20 mil in series with the model inductors. 

5.3 Design of Model Self Inductors 

The inductance of a toroid is given by 
4N^ M A X 10"^ 

L = g henries, (5-10) 

where 



N = number of turns 
)j^= relative permeability 

2 

A = core cross sectional area (cm) 

D = mean core diameter (cm) . 

For the Arnold core A-21 8259-2 used in this design, = 205, so that 
L = 259 mh per 1000 turns on any full wound core. A full wound core 
is defined as one in which the winding occupies three-fourths of the 
available window area with a packing factor of 0.7. For any other 
full wound coil on the A-21 8259-2, 



N = 1000 




(5-11) 



where "L" is the desired inductance in mh. The inductances shown in 
Table 4.5 require the turns shown in Table 5.1. The listing of turns 
required in the 8-mile sections is shown in Table 5.2. To achieve 
the Q characteristic shown in Figure 5.1, it is necessary to have the 
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Table 5.1 

Turns Required in 4-Mile t.l.m. Inductors 









93 


91 


90 


94 


92 


91 


95 


93 


92 


96 


94 


93 


97 


95 


94 


98 


96 


95 




Table 5.2 




Turns 


Required in 8-Mile t.l.m. 


Inductors 








131 


129 


127 


133 


130 


129 


134 


131 


130 


136 


133 


131 


137 


134 


133 


139 


136 


134 



core full wound. In the case of the 4-mile inductors this works out 
to be (No. 10 AWG, and for the 8-mile sections - No. 11 AWG. For the 
4-mile sections, a nominal twenty feet of No. 10 is required at 1.0 
mfi/ft for 20 mJ^ per inductor. For the 8 mile sections, a nominal 
twenty-eight feet of No. 11 is required at 1.26 mfi/ft for 35 mfi 
per 8-mile inductor. In Tables 5.1 and 5.2, the first five numbers 
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in each column are the turn numbers at which a tap must be provided. 

Referring to equation (5-10) it is seen that the inductance 
of a toroid is independent of current if is independent of 
current. A graph of vs. a-c flux density taken from reference 18 
is reproduced in Figure 5.3. It is seen that p^ remains within five 
percent of nominal up to 4000 gauss. From basic inductor theory: 



= i = NBA_ . NpHA 



" I " I " I " I 



(5-12) 




(5-13) 

(5-14) 

(5-15) 



ttDH = NI 



H = NI/ttD 

Substitute (5-15) into (5-12): 



NBA ^ N^pA 
I ttD 



(5-16) 



i = M 

I nD 



(5-17) 



ttDB 
■ Np 



(5-18) 



In equation (5-18) 



I = complex amplitude of current (amperes) 

N = no. of turns in winding 
D = mean diameter of core in meters = .0463 m 
B = peak flux density (wb/m ) 
p = normal permeability = Pj,Pq* 



For the t.l.m. inductors at 4000 gauss: 



B = .4 wb/m^ 

p = 205(4 ttx10'^)^ = 2.58x10 
D = 0.0463 m 



N = 93 or 133 . 
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So that I peak = amperes 
=1.70 amperes, 

and the inductance of the t.l.m. inductors remains within five percent 
of nominal up to 1.70 amperes r.m.s., or 31 p.u. 

5.4 Design of Model Mutual Inductors 

Figure 2.2 shows the mutual inductors and L 23 (L|^^) 

as transformers between the phases. Since the mutual inductances are 
small, and since there is sufficient space remaining in the windows 
of the self inductors, is wound over L^, and is wound over Lj^. 
The mutual inductances to be provided in the 4-mile sections are 
shown in Table 4.6; those for the 8-mile sections are simply double 
the values in Table 4.6. From basic mutual inductor theory (See e.g. 
reference 15 p.p. 178 ff) it is knov/n that, for example. 



'■ab " *^\/^^ab^s '"a 



and from equation (5-11) for the model core. 
Where (all in mh) : 



(5-19) 

(5-20) 



(5-21) 



Lab ~ mutual inductance between phases "a" and "b 
(L^l^)g= self inductance of 

k = coefficient of coupling between and 
L" = measured mutual inductance when and 
are connected to reinforce. 

L" = measured mutual inductance when and 



II 



are connected to oppose. 
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(Nab)$ = of turns to obtain (L^b^s* *"ab^' 

Several tests were performed to determine k by winding six, seven, and 
eight turns of (L^b^s ^a ^ eductor of 2.49 mh. Knowing the 

number of turns, (L^b^s obtained from (5-21), and L^b was obtained 
as indicated by (5-20). (Lgb)s ®od L^b were then substituted into 
(5-19) to determine k. Under these conditions, k for the t.l.m. 
inductors was found to be 0.986. Solving (5-19) for (Lab)s ®od 
substituting into (5-21), 



<"ab>s = « 




(5-21) 



With taken as 2.49 mh, 

a 



<Nab>s ’ 39.8 (5-22) 

and the number of turns to achieve the mutual inductance values, L^b* 
of Table 4.6 are 6, 7, 8, 9, 10, and 13, so that taps must be provided 
at 6, 7, 8, 9, and 10 turns. With the taps set, the self inductance 
corresponding to each tap can be found using equation (5-21) and then 
the values of L^b other values of can be found using (5-19); 
the results are shown in Table 5.3. A similar procedure applied to 
yields. 



<\c>s = “O-S l-bc 

and still the taps are at 6, 7, 8, 9, and 10 turns, with the 13th 
turn providing Lb^ = 0.32 mh. A list of Lb^ values for a 4-mile 
section is given in Table 5.4. In order to wind L^b ^bc 
Lg and Lb respectively, it is necessary to use four feet of No. 12 
AWG magnet wire at 1.59 mfl/ft for an additional resistance of 6.4 mQ 
in series with the phase "b" and "c" lines. A list of mutual 
inductance values to be built into the 8-mile sections is shown in 
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Table 5.5. 



Table 5.5 



ab'"*’) 


Lbc(inh) 


0.30 


0.30 


0.36 


0.36 


0.40 


0.40 


0.46 


0.46 


0.50 


0.50 


0.66 


0.56 



Table of Mutual Inductance Values 
Designed into the 8-Mile Sections 

The taps to yield (using = 4.98) are 8, 10, 11, 13, 14, and 

18; taps for (L| 3 j .)3 (using = 4.82) are 8, 10, 11, 13, 14, and 17. 
The corresponding values of (Lgij)^ for the 8-mile sections are (in mh) 
0.0166, 0.0259, 0.0312, 0.0435, 0.0505, and 0.0834. For 
0.0166, 0.0259, 0.0312, 0.0435, 0.0505, and 0.0745. These tabulations 
along with those of Tables 5.3 and 5.4 show that the self inductance 
of a mutual winding is negligible when compared to the self inductance 
of the line. In addition. Tables 5.3 and 5.4 indicate that the change 
in or resulting from a change in or L|^ over their entire 
range is less than five percent and therefore negligible also. In 
order to wind and over and in the 8-mile sections, it 
is necessary to use approximately six feet of No. 12 AWG magnet wire 
at 1.59 mfl/ft for an additional resistance of 9.5 mft in series with 
the phase "b" and "c" lines. 
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CHAPTER 6 

CAPACITOR SELECTION 

6.1 Factors Determining Capacitor Selection 

As discussed in detail in Chapter 4, eight diverse transmission 
line types were analyzed in the EPSEL Line Parameters Program, and 
these eight types were taken to cover the range of Typical EHV lines. 
The model capacitance values implied by these lines for a 4-mile model 
section are shown in Tables 4.4a and 4.4b. A complete listing of the 
capacitance parameters to be covered by both the 4 and 8-mile model 
sections is shown in Table 6.1. Having scaled the eight representative 
lines to result in a minimum inductance range, the range of each 
capacitance is seen to be somewhat larger. It is proposed to cover 
the range of each parameter of Table 6.1 with a bank of six capacitors 
wired and brought to the operating panel so that any combination of the 
six capacitors in each bank can be arranged in parallel. Thus, with 
one or two capacitors to set the lowest value in a range, approximately 
five or four capacitors are available to cover the remainder of the 
range. Thus there are a maximum of 120 and minimum of 24 capacitance 
values available to cover the range of each capacitance parameter 
indicated in Table 6.1. With this "dense" provision of model values 
it is possible to set any value in the span within a very few percent 
('^3%) if the tolerance of the capacitors is less than a few percent. 

It is proposed to provide a two to three percent tolerance on any 
capacitance selected by the operator so that the over-all capacitance 
provision of the model is within about 5% of the value required to 
model any given line. With expected transient overvoltage maxima 
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of about four p.u., all capacitors must be capable of withstanding short 
duration peaks up to 500 or 600 volts. In addition to tolerance and 
voltage requirements the capacitors must have extremely low losses. 

An indication of tiie dissipation factor required of the model 
capacitors can be obtained by a simple analysis of the "dissipation 
factor" of a typical line. Considering the positive sequence parameters 
of a type 41A line using the values given in Table 4.1, 

S 

Cp . = = .0181 yf/mi. 

Cp 2 = .0725 yf/4 mi. 



Rp ^ = .216 ft/4 mi. 

Using the "series model" for a lossy capacitor: 

Op 2 /4 mi - wRp ^ Cp ^ (6~1 ) 

Dp 5 = 0.0000059/4 mi. 0 60 Hz (6-2) 

Dp 3 = 0.0000985/4 mi. 0 1000 Hz (6-3) 

That the values given by equations (6-2) and (6-3) are the values 
required of a 4-mile model section can be seen by a corresponding 
calculation for the model. Using Table 4.4a for the type 41A model: 

T T C + Cl + C 

'^P.S. - S * — T * ‘•ab * <^30 * Sc 



Cp ^ = .230 yf/4 mi 
R 



j R^ + R. + R^ 

p 5 4 ^ — 3 ^ ^ = -068 ft/4mi 



Dp 5 = 0.0000059/4 mi . O 60 Hz 

Dp 3 = 0.0000985/4 mi. P 1000 Hz. 



(6-4) 

(6-5) 



p? 
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at low currents. The basic relation governing the Q of an inductor can 
be written as (for the frequencies of interest here) 



Q = 



o)LI‘ 



R + R..I^ ’ 
dc ac 



Where 



(5-1) 



L = Inductance in henries 
w = 2n X frequency 

Rdc “ dc resistance of the winding (n) 

R._ = equivalent core loss resistance (fl) 

dC 

2 

With these definitions it is seen that R_^I corresponds to both the 

dC — — 

hysteresis and eddy current losses (See reference 17 p.p. 58ff.) and 
thus represents the total core loss. R,_, the core loss coefficient, 

d C 

1 8 

is itself a function of current. For powdered cores the manufacturer 
provides data on R explicitly, whereas for Fe-Si the entire term 

d C 

2 

R .1 is provided. If the data provided by the manufacturer is 
dC 

converted to the terms of equation (5-1), it is found that for Fe-Si 
the core loss, W, is given by (approximately) 

W = k^(w)I^•^ (5-2) 

whereas for the powdered Mo-Ni the loss is given by (approximately) 

Vi = Rgj.1^ = l<2((o)I^-^ . (5-3) 

In equations (5-2) and (5-3), k^(w) and k 2 (w) are core loss constants 
and are a function of frequency. If equations (5-2) and (5-3) are 
substituted in equation (5-1) there results, respectively. 



Q 



Fe-Si 



ojLI^ 



(5-4) 



and 



M 
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occurs primarily at low frequencies in connection with switching opera- 
tions which connect the terminals of a charged line to large inductive 
reactances. Natural frequencies of the order of 40 Hz are known to occur 
under these conditions with over-all Q approaching 300. 

Host capacitor manufacturers specify or test dissipation factor at 
1000 Hz, so that at 1000 Hz it is possible to measure capacitance and 
compute resistance using either the parallel or series model of a lossy 
capaci tor: 



D = 



_! 



( 6 - 6 ) 



The values of or Rp so obtained, however, are true only for a limited 

range of frequencies and cannot be used at 40 Hz. As frequency decreases 

toward zero the value of Rp (resistance in the parallel model) approaches 

the insulation resistance of the dielectric - which for polystyrene is 

20 

of the order of ten million megohms for a 0.15 pf capacitor. The 

dissipation factor, D, of a typical .15 yf polystyrene capacitor at 

20 

1000 Hz is given as 0.0001. This value of D implies Rp = 10.6 megohms 

at 1000 Hz. Few materials other than Teflon or a modified polystyrene 

can equal this low value of D or high value of Rp at 1000 Hz. If 

equation (6-6) is solved for Rp in terms of w, R^, Cp, and 

(C^ = Cp), and set equal to 1.0 mn, = 0.15 pf, and f = 40 Hz the 

value of Rp obtained is 4.4 million megohms. Thus in order to obtain 

R^ as low as .001 n, it is necessary to have Rp as high as 4.4 million 

megohms at 40 Hz. Little data exists on any dielectric material below 

21 

100 Hz - standard reference works supply no data below 100 Hz on any 
material and give no information vs. frequency on polystyrene. In an 
effort to obtain an estimate of Rp for the 0.15 pf polystyrene capacitor 
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at 40 Hz, it was assumed that Rp remains constant at the 10.6 megohm 
value down to 100 Hz and then rises linearly to 10 million megohms at 
d.c. This procedure gives Rp as 6 million megohms at 40 Hz. (This is 
a conservative estimate in the opinion of the author.) No material 
less costly than polystyrene can approach within an order of magnitude 
of this value. Polystyrene capacitors are therefore selected to form 
the "base" capacitance of each capacitor bank in the model. 

6.2 Capacitors Selected for Model Sections 

Since the polystyrene capacitors conservativety (but not by 
orders of magnitude) meet all line requirements, a compromise in favor 
of lower cost is made in selecting the trimming capacitors in each 
bank. Thus, while polystyrene capacitors will be made to provide 
80 percent (approximately) of any capacitor value selectable by the 
operator, the remaining 20 percent (approximately) is to be made up of 
lower cost mylar or silvered mica capacitors. The latter capacitors 
have losses ten to fifty times higher than the polystyrene but they 
constitute only ten to twenty percent of a selected capacitor value 
and so should not appreciably lower the high quality provided by the 
polystyrene base capacitors. The severe low-loss constraint applies 
mainly to the 765 kV line, so that the polystyrene capacitors in each 
bank are selected so that for the 765 kV case they will constitute no 
less than 90% of the required capacitance of a 4-mile section. The 
polystyrene capacitors selected for the t.l.m. are Arco Manufacturing 
Co., type PJ, 400 VDC, +2%, tubular capacitors. The 400 V units 
were judged adequate since each capacitor receives a minimum of two 
sixteen hour burn-in periods with 300 percent of rated voltage applied 
at 85®C. The mylar capacitors are Cornell Dubilier Co., type WMF6, 
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600 VDC, + 10 percent, tubular capacitors. Two small capacitance values 
below the range of standard mylar sizes (330 pf and 680 pf) are Cornell 
Dubilier Co., type CD19FD, silvered mica, 500 VDC, ^ 5%. The capacitor 
types and sizes selected to span each capacitance parameter (the 
capacitor banks) are shown in Table 6.2. 
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CHAPTER 7 

RESISTANCE AND LOSS CONSIDERATIONS 
7.1 Over-all Q Considerations 

As discussed in Chapters 4 and 5, the Q requi rements of this 
t.l.m. are dictated by the need to model the 765 kV line. This con- 
straint is particularly severe in the 4-mile model sections since 
connection resistance is expected to take up any margin provided by 
designing the inductors with double the required Q at 60 Hz. In 
Chapter 4 it is indicated that the positive sequence Q of a typical 
765 kV line is twenty. The resistance required in each phase to achieve 
this result is shown in Table 4.4a - 37 mn in phases "a" and "c" and 
43 mfi in phase "b". As discussed in Chapter 5, 20 mfi is "provided" in 
each phase by the resistance of the wire to wind the self inductors 
La* L|^, and L^. An additional 6 mfl is provided in phases "b" and "c" 
due to the fact that these phases contain the resistance of the windings 
for and respectively. (Actually, and for the typical 
765 kV line are provided by the 7-turn-tap so that the resistance 
added by the and L|^^’ windings will be approximately 4 mil in the 
4-mile sections. The operating panel resistance to be provided at this 
point in the accounting is 17 mil in phase "a", 19 mil in phase "b", and 
13 mil in phase "c". But connection resistance has not yet been 
accounted for in these figures. 

In a typical model set up (See Chapter 8), six to eight banana 
jack connections will be made in series in each section from input to 
output, so that the resistance of these must be accounted for. The 
jacks selected for this design are E.F. Johnson Co., standard banana 




0 
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jack, cadmium plated insert, catalog no. 108-0900-001. The corresponding 
plug is E.F. Johnson Co., standard plug, catalog no. 108-0300-001. Six 
of these jacks and plugs were connected in series as shown in Figure 7.1. 
Sixteen inches of solid, no. 14 AWG, copper bus bar wire was used, and 
all connections were heavily soldered. The total resistance was measured 
at 12 mJ2 with the GR 1650A resistance bridge. The resistance of the 
no. 14 wire is 3.3 m^2 so that the resistance of each connection averages 
to about 1.5 mJ2. All plugs and jacks were new, clean, tight fitting, 
and all connections were made and remade several times - always with the 
same result. The resistance accounting to this point is: 

Rg = 8 mft, = 5 mfl, R^ = 1 mil. In other words, at this point the 

resistance values required at the front panel to bring the section 
positive sequence Q down to twenty at 60 Hz are the values cited just 
above. Each 4-mile section; exclusive of external hookup wire, but 
including operator panel connector resistance is seen to have built 
into it. 



To provide an exact positive sequence Q of twenty at 60 Hz, it is 
necessary to provide for operator selection; R^ = 8 mil, R^j = 5 mfl, 
and R^ = 1 mQ. As pointed out in Chapter 5 in connection with 
Figure 5.1 however; setting the Q to twenty at 60 Hz for a 765 kV line 
model will result in the model Q falling below line Q above about 
1700 Hz. Thus, it may be desirable under certain circumstances to set 



Rg (intrinsic)^ = 29 mn 
% 

R|^( intrinsic)^ = 38 mn 

R^(intrinsic)^ = 38 mil 
c 4 



(7-1) 

(7-2) 

(7-3) 
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Test Circuit for Measuring Resistance 
of Banana Jacks and Plugs 



Figure 7.1 
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Q at 60 Hz slightly above twenty to ensure that the model Q remains 
above line Q up to 3060 Hz. 

An accounting procedure identical to that above applied to 
the 8-mile model sections gives the results. 



The values of R,, R. , and R^ required in the 8-mile sections to model a 

a D C 

765 kV line are (doubling the appropriate values in Table 4.4a): 



It is seen that the 8-mile sections have considerably more margin than 
the 4-mile sections and, in fact, using the intrinsic values of 
equations (7-4)-(7-6), the positive sequence Q of the 8-mile sections 
can be set as high as 31 at 60 Hz. 

7.2 Resistance Selection 

Referring to Table 4.2 (excluding the 1500 kV lines) the values 
of positive sequence Q to be provided are seen to range from about 
eight to twenty. To achieve the Q values lower than twenty, resistance 
values must be provided at the operating panel for insertion with the 
intrinsic resistance of the section. For the 4-mile sections, the 
lowest resistance value to be made available is clearly zero since 
this is the required value to ensure positive-sequence Q's just over 
tv/enty. For the 8-mile sections a similar zero base will be provided 
since it would simply not make sense to deny the operator the selection 




(7-4) 

(7-5) 

(7-6) 




(7-7) 

(7-8) 

(7-9) 
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of positive sequences Q's to thirty should the occasion arise. A list- 
ing of Q values to be provided for each phase is shov/n in Table 7-1. 



Table 7.1 

4-Mile Section 8-Mile Section 



Qa 


Qb 


Qc 


Qa 


O' 


Qc 


30 


22 


21 


40 


30 


28 


25 


19 


18 


34 


25 


24 


20 


16 


15 


28 


20 


20 


15 


13 


12 


22 


16 


16 


12 


11 


10 


16 


12 


11 


10 


9 


8 


11 


9 


8 



Table of Q Values Provided 
in the T.L.M. Sections 



The operating panel resistance values to achieve the Q's shown in 
Table 7.1 are shown in Table 7.2. 



Table 7.2 

4-Mile Sections 8-Mile Sections 



R 


K 


R 


R 


K 


R 


a 


b 


c 


a 


b 


c 


0 


0 


0 


0 


0 


0 


6 


7 


8 


8 


11 


12 


15 


15 


17 


19 


28 


25 


30 


27 


30 


36 


49 


46 


44 


40 


43 


66 


85 


92 


59 


67 


64 


116 


132 


147 



Table of Operating Panel Resistance Values 
Corresponding to Q's of Table 7.1 (R in mJ2) 
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The resistance values indicated in Table 7.2 were obtained using 

o 

Ri - R^- (intrinsic) (7-10) 

(Lg)^ = 2.34 mh, (^)^ = 2.26 mh, (L^)^ = 2.18 mh 
(Lg)g = 4.68 mh, (Ljj)g = 4.52 mh, (Lj,)g = 4.36 mh 
i = a,b,c for each section length 
Q. is the corresponding Q value from Table 7.1. 

vJ 

It is seen in practice, rather than setting positive sequence Q per se, 
the Q of each phase is set individually using Tables 4.4, 7.1, and 7,2. 
The resistance indicated in Table 7.2 is achieved by connecting 
resistances in series so that each indicated value is the sum of the 
resistances required to form each of the preceeding resistors. The 
actual resistors required to obtain the operating panel values indicated 
in Table 7.2 are shown in Table 7.3. 

Table 7.3 



4- 


•Hile Sections 


8-1 


'lile Sections 


, 's 


Rk's 


R 's 


R 's 


Rk's 


R '! 


a 


b 


c 


a 


b 


c 


6 


7 


8 


8 


11 


12 


9 


8 


9 


11 


17 


13 


15 


12 


13 


17 


21 


21 


14 


13 


13 


30 


36 


46 


15 


27 


21 


50 


47 


55 


: of 


Resistors 


in 


the T.L 


.M. Sections 



The resistors indicated in Table 7.3 are to be "manufactured" by 
cutting appropriate lengths of nickel or copper-nickel alloy ("cupron") 
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resistance wire. A very large selection of resistance wire sizes and 
types are available from the Driver Harris Co. of Newark, a type 

particularly suited to the smaller resistors above is the company's 
0.035" diameter nickel wire of 0.045 ohms per foot. For the larger 
resistors no. 20 AWG "Cupron" (copper-nickel) wire with a nominal 
resistance of 0.287 ohms per foot can be used. The latter wire can 
safely carry a direct current of three amperes indefinitely but is 
hot to touch at this current level. Three rms amperes is far in 
excess of that required to simulate the most severe faults, so that 
no. 20 AWG resistance wire should be adequate for this t.l.m. 




1 
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CHAPTER 8 

MODEL PHYSICAL CONSTRUCTION 

8.1 Electrical Construction of a Typical Section 

As discussed in Chapter 4, this t.l.m. consists of thirty 4-mile 
sections and thirty-five eight mile sections. A simplified schematic 
diagram of the wiring of a typical section is shown in Figure 8.1 and 
is the same for both the 4-mile and 8-mile sections. The component 
values are not indicated since they are different for the two section 
types, but the required values can be taken from the tables of Chapters 
5, 6, and 7. In wiring each section, all hookup wire in a direct path 
from input to output is No. 10-stranded, tinned, copper wire. Capacitor 
hookup is effected with No. 16 stranded copper wire. All connections 
are soldered. 

8.2 T.L.M. Physical Construction 

Each t.l.m. section is constructed on the chassis shown in 
Figure 8.2. The sixty-five chassis are made of plexiglass of the sizes 
indicated - the operator's panel is of black plexiglass. The layout 
of the operator's panel is shown in Figure 8.3. The jacks are color 
coded by "phase" and the circuit diagram is painted by a silkscreen 
stenciling process. Each of the sixty five chassis is assembled in a 
larger console as shown in Figure 8.4. The chassis slide into the 
console on rails and are secured with screws through the operating 
panel. The chassis are inserted into the console from both "front" 
and "rear", so that twenty-four sections are accommodated per console. 
Three consoles such as that shown in Figure 8.4 are required with a 
fourth "half-console" for the output section and terminal equipment. 
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T.L.H. Console 



41" 



Section 
No. 1 


No. 12 
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No. 3 
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CHAPTER 9 

SUMMARY AND RECOMMENDATIONS 

9 . 1 Summa ry 

This report presents in detail the considerations in the design 
of a transmission line model for use in the M.I.T. Transmission System 
Simulator. The design includes the modeling of the inductive effects of 
non-transposi tion by incorporating two linear mutual inductors between 
phases. The inter-phase resistance effects are approximated by assuming 
that the three inter-phase resistances are equal. A circuit diagram of 
one three-phase L-section is shown in Figures 2.2 and 8.1. Note that 
the mutual inductors are physically wound on the respective self inductor 
core and that the return circuit is provisional pending a complete 
design due in September, 1971. 

In order to deduce the number of sections required for the model 
and establish a quantitative measure for predicting the performance of 
model types with various terminations, the pole error analysis of Chapter 
3 was developed. This analysis is based upon a comparison of the poles 
of the driving point impedance of actual lines with those of proposed 
models. Solutions for the poles of tt and L-section models were obtained 
by computer for models up to forty sections with open and short circuited 
terminations. The pole error curves of Chapter 3 show in quantitative 
terms the convergence of the poles of various models toward those of 
continuous lines with similar terminations. A technique is outlined 
(See Figure 3.11) in which the pole error curves are used to deduce 
the number of sections required to meet any arbitrary accuracy spec- 
ification. The specification set for the M.I.T. model is that the 
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error in the pole frequency of the model shall be within five percent 
of the frequency of the corresponding actual pole over the range from 
sixty to 3060 Hz for line lengths of sixty to four hundred miles for 
both open and short circuited terminations. This criteria led to the 
decision to construct thirty four-mile model sections and thirty-five 
eight-mile sections ("L" configuration). 

In Chapter 4 parameter data for eight representative EHV line 
configuration types are presented in order to deduce the parameter values 
required in the four and eight-mile model sections. Several scaling 
techniques are suggested for scaling the parameters of the diverse line 
types to the parameters provided in the model. Careful scaling suggests 
that only a narrow range of component values need be provided for each 
parameter in the model. In the M.I.T. design, six values are provided 
for each of the sixteen parameters to be modeled. 

The quality factor (Q) of typical EHV lines are seen (Figure 4.1) 
to vary markedly with frequency in the range from sixty to 3060 Hz. 
Positive sequence Q increases monotonically to over 175 for the highest 
Q lines studied, and zero-sequence Q maximizes near 1000 Hz and 
generally remains constant up to 3060 Hz. The requirement to model 
these Q characteristics is a major factor in the design of all model 
inductors and in the selection of all model capacitors. The positive 
sequence Q characteristics is achieved by the use of molybdenum- 
nickel toroids for all phase inductors. In fact, the model inductors 
have an intrinsic Q characteristic approximately double the value 
required to model the highest Q EHV lines in operation. Preliminary 
plans for modeling the return path or zero-sequence inductance and 
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Q suggest a separate modeling of the ground wire path and earth return 
path - the final design is to be completed by September, 1971. 

A careful account of wire and connector resistance shows that 
for the section lengths chosen, every milliohm significantly reduces 
the intrinsic Q provided by the model inductors. To avoid excessive 
losses and provide for operator adjustment of Q, all connections are 
soldered, and all wire on a direct path from input to output is 
no. 10 AW6 hook up wire. To avoid further losses, and to model the 
condition of a charged line isolated with a high Q reactor at its 
terminals, extremely low loss polystyrene capacitors are incorporated 
into each t.l.m. capacitor bank. Approximately eighty-five percent 
of any model capacitance parameter selectable by the operator is low 
loss polystyrene in parallel with approximately fifteen percent 
additional "trimming" capacity provided by mylar units. Overall 
capacitance tolerance in the t.l.m. is within five percent of any value 
selected. Six capacitors are provided for each capacitance parameter 
and all capacitors are available for operator connection into any 
desired configuration. 

The final result for the M.I.T. t.l.m. is a total of sixty-five 
L-sections for modeling line lengths to 400 miles. Average Q in the 
four-mile model sections is a maximum of 24 at 60 Hz; the maximum 
average in the eight mile sections is 32. All capacitance and 
inductance selection within the range of a parameter is within five 
percent of any desired value. A similar five percent tolerance is 
planned for the return path inductance and Q. Q values less than the 
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maxima given above are obtained by the insertion of appropriate resistors 
provided in each section. 

9.2 Recommendations 

The design of this t.l.m. has progressed in parallel with its 
actual construction and component procurement. This was attempted in 
view of a master's thesis time constraint. The basic scale factor 
between the model and lines was selected in order to make the t.l.m. 
compatible with existing laboratory equipment. In fact, the t.l.m. 
ultimately will represent an investment several times the cost of the 
TSS transformer which dictated the scale ratio. This scale ratio 
(the 0.355 for the 765 kV line) was selected and found acceptable 
primarily on the basis that it resulted in a feasible, low cost inductor 
core without scaling the capacitors out of the range of economical 
values provided "off the shelf" by manufacturers. Not until after the 
cores had been purchased was the significance of capacitor dissipation 
factor fully appreciated. Although the required capacitor values are 
well within the range of normal capacitors, they are among the highest 
and most costly polystyrene units. The result was the requirement to 
purchase capacitors nearly three times the cost of inductor cores. 
Polystyrene capacitors are required regardless of what reasonable scale 
ratio is selected, but it is suggested that perhaps a better cost trade 
off could have been achieved had the TSS transformers been redesigned to 
permit a scale ratio closer to unity. This would require larger, more 
expensive, cores for the t.l.m. inductors but would significantly reduce 
the size of capacitors required. The larger polystyrene units are very 
expensive, so that even a small reduction in values required could result 



I 




118 



in considerable savings. This savings could be directed toward purchasing 
a better (larger) inductor core permitting even more margin in Q by making 
the connector resistance a smaller fraction of the total resistance of the 
section. It is recommended that, should a similar project be undertaken 
in the future, the entire design should be completed in every detail, a 
prototype section constructed, and a cost analysis completed before a 
single large scale purchase is made. With a better cost trade off between 
t.l.m. components and a more complete view of overall t.l.m. investment 
compared to other TSS components, a more informed decision can be made as 
to what part of the TSS should be scaled to achieve compatibility between 
them al 1 . 

Concerning the pole error analysis of Chapter 3, several comments 
seem appropriate. It should be noted that none of the pole error analysis 
was carried to more than 39 sections. This is due to the fact that 
none of the root subroutines available was able to generate stable 
solutions beyond this point. Several algebraic root finding methods were 
used in an effort to obtain solutions above the 40th degree, but none 
succeeded. This is believed due to the fact that all the roots are in 
the region between minus one and plus one, and as the degree of the 
polynomial gets higher the roots tend to crowd near the extremes. Scal- 
ing of the coefficients of the polynomials was attempted but this 
technique resulted in instabilities at even lower degree. Three of the 
four subroutines used in this problem are available in the IBH System/360 
Scientific Subroutine Package (360-CM-03X version III) using three 
different (but algebraic) techniques. The subroutine POLRT uses the 
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Newton- Raphs on iterative technique and is explicitly limited to poly- 
nomials of degree thirty-six or less. The subroutine DPRQD uses the 
quotient-difference-displacement method - this routine generally 
became unstable at about the 30th degree for this problem. The DPRBM 
subroutine uses the Bairs tow iterative method of quadratic factoriza- 
tion; this routine was the most stable of the four methods - giving 
stable solutions through degree thirty-nine. A fourth subroutine, 

ROOT 1, using Muller's method (written by the computation group at 
M.I.T.) became unstable at about degree twenty-four. The Bairstow 
method is the routine selected for all the programs of Appendix D. 

No extensive effort was made to determine exactly why the subroutines 
became unstable, nor were more elegant matrix methods attempted. It 
was felt that solutions to degree forty gave information sufficient 
for the purposes of this design, so that no further effort was expended 
in this direction. Nevertheless the instability problem is at least 
of academic value and worth further study. 
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APPENDIX A 

MODEL TRANSFORMATIONS 



A.l Equations of the General Non-Transposed Model 

In the general, non- transposed model of Figure 2.1 Chapter 2, 
the differential equations relating the line to ground voltages and 
the currents are given by Kirchoff's Voltage Law (The Voltage Law): 



3V. 

C 

dz 



3z 



9Z 



31 . 



3L 



^^'ll'^^OO^^a ^^ll‘^‘-00^ 3t ^^OO^b ^^12‘^'-00^ 3t 

31. 



(L-iQ+Lp,/^) 



00"c '13 00' at 



31 



31. 



%0^a ^ ^'- 12 ‘^'- 00 ^ at ■^^‘^22'^'^OO^^b ^‘-22'^'-00^ at 

31 



31. 



*^00^c ^‘-23'^‘-oo^ at 

31, 



^OO^a * ^‘-13‘^‘-00^ at^ * *^ 00^0 ^*-23‘^'-00^ at" 



(A-la) 



(A-lb) 



ai 



^'^33'‘‘^00^^c ^*-33'‘''-00^ TT 

Similarly, the currents can be related to the line to ground voltages by 
Kirchoff's Current Law (The Current Law): 
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A. 2 Transformation of the General Non-Transposed Model to a "Balanced 
Mutual" Model 

Equations (A-1 ) describing Figure 2.1 can be transformed to a 
set describing Figure 2.2 using simple algebra. Let ~ 3 Iq 
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Equations (A-2) and (A-3) are the relations describing the "balanced 
mutual" model of Figure 2.2 with: 
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A, 3 Transformation of the General Non-Transposed Model to a Transposed 
Hodel 



One transposition cycle of a transposed three-phase transmission 



line can be visualized physically as shown in Figure A.l. In the 
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Figure A.l 
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transposed line let: 

lL = self inductance of a phase = 

lTj = mutual inductance between two phases = 

T T 

= capacitance to ground of a phase = Cq 

T T 

= mutual capacitance between two phases = 

r|^- = self resistance of a phase = 

Then from figure (A-1 ) it is seen that in terms of the parameters of 
Figure 2.1 : 



+ Loo + L, 



11 


‘■22 ^ 


'-33 


1 1 


cc 

3 




v)0 J 1 

■-s 


(A-5a) 


T 




'-L " 


'-12 


*h3 


+ 


'-23 , T 


(A-5b) 


12 ■ 




3 




“ ‘-M 


T 


CT . 
^20 


O 

CO — 1 

o 

II 


^10 


^20 


+ 


^30 pT 


(A-5c) 


10 ■ 




3 






T 


cT = 
h3 


cT = 
^23 




+ C,3 


+ 


^23 pT 


(A-5d) 


12 ■ 




3 




" ni 


J _ 


= 

*'22 


rT - 

^^33 


'^11 


+R 22 


+ 


^^33 _ .T 


(A-5e) 


11 ■ 




3“ 







The significance of these relations can be seen more clearly by modeling 
a transposition cycle of the transposed line with the model of Figure 
2.1. The differential equations corresponding to equations (A-1) with 



I + I. + I = 3K are: 
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Adding and subtracting LjJj — as appropriate in equations A -6 obtain: 
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In the balanced steady-state Ig = 0; the currents in the transposed line 
are positive sequence only, and equations (A-7) reduce still further to 
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three decoupled "single-phase equivalent" equations. With L^-LjI|=Lp ^ = 
positive-sequence inductance and ^5 = Rp 5 = positive-sequence resistance 
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If only zero-sequence currents flow in the line, I^ = Ij. 
equations (A-7) reduce to 
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Equations (A- 8 ) and (A-9) describe (in terms of The Voltage Law) the 
symmetrical components equivalent circuits of the transposed line 
described by equations (A-7). The transposed line corresponding to 
equations (A-7) is shown in Figure 2.3. 
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The Current Law equations for the transposed line corresponding 
to equations (A-2) can be seen by inspection of Figures 2.1 and A-1 : 
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In the balanced steady-state Vq = 0; the voltages are positive sequence 
only, and equations (A-11) reduce to three decoupled "single-phase 
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If only zero-sequence voltages are applied to the transposed line, 
Va = Vb = Vc* And with S ~ sequence capacitance, 

equations (A-11) become: 
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Equations (A-12) and (A-13) 


describe (in 


terms of The Current Law) the 



symmetrical components equivalent circuits of the transposed line 
described by equations (A-11). The transposed line model corresponding 
to equations (A-11) is shown in Figure 2.3. 
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APPENDIX B 

MODEL IMPEDANCE DERIVATIONS 

B.1.1 Derivation of L-Section Model Impedance- Inductor Input- 
Shorted Terminatfon 

UsinQ ths ni ssction lin6 shov/n in Figur6 3.3, writs diffsrencs 
equations near the nth section (0<n<m): 
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Subtracting (B-1 ) from (B-3): 
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Next assume e (where X may be a complex constant, 

X=A^+jA^. ) and substitute into equation (B-6): 
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Therefore, provided X satisfies (B-7), the solutions to equation (B-5) 



take either of the forms: 

Vn = Re[V|^ e^^] = Re[(K+ 0 “^^+ K. e'"^)e^^] (B-8) 

= Re[(A cosh nX + B sinh nX)e^^] (B-9) 

To obtain the current, I^, use equation (B-1) assuming solutions of the 



form: 
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In terms of equation (B-9) equation (B-12) can be written as: 

T _ A(cosh nX - cosh(n-l)X) + B(sinh nX - sinh(n-l)X) 
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If the short is connected at the zeroth node, then Vq = 0 and from 
equation (B-9) it is seen that A = 0, so that elsewhere on the line: 



V = B sinh nX . 
n 

And from equation (B-13): 
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Finally, the impedance looking into the mth node 
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Where it is noted again that X must satisfy 
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(B-7) 



In the case that the model is lossless (R^ = 0), equations (B-16) and 
(B-7) reduce to: 



The solution for the poles of (B-17) is indicated in Appendices B.1.2, 
C.2.1 , and D.l . 

B.1.2 Solution for the Poles of the L-Section Model-In ductor Input- 
Shorted Terminati^ ~~ — 

The solution for the poles of equation (B-16) or (B-17) rests 

upon recognition of the fact that the denominator of (B-16) and (B-17) 

is a polynomial in cosh X of degree (m-1) for all m. That this is 

true can be seen by expansion of the denominator for the first few 

values of m. The expansions of sinh mX up to m = 7 are shown in 

Table B.l . 
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Table B.l 



Expansions of sinh mX 



sinh X = sinh X 



sinh 2X = sinh X(2 cosh X) 

sinh 3X = sinh X(4 cosh^X-1 ) 

sinh 4X = sinh X(8 cosh\-4 cosh X) 

sinh 5X = sinh X(16 cosh^X-12 cosh^X+1 ) 

sinh 6X = sinh X(32 cosh^X-32 cosh^X+6 cosh X) 

sinh 7X = sinhX(64 cosh^X-80 cosh^X+24 cosh^X-1) 



131 



The expansions listed in the table were obtained by repeated 
application of the identity sinh mX=sinh(m-l )X cosh X+cosh(tn-l )X sinh X. 
The use of this identity to obtain sinh tnX requires the preceeding 
recursion terms sinh(m-l)X and cosh(m-l)X. The term sinh(m-l)X is 
simply the preceeding term in the table, but cosh(m-l)X must be obtained 
from a separate process of expansions of cosh mX. Table B.2 shows the 
first few expansions of cosh mX generated in the loop required to get 
sinh mX. 

Table B.2 

Expansions of cosh mX 
cosh X = cosh X 
cosh 2X = 2 cosh^X-1 

3 

cosh 3X = 4 cosh X-3 cosh X 

cosh 4X = 8 cosh^X-8 cosh^X+1 

cosh 5X = 16 cosh^X-20 cosh^X+5 cosh X 

cosh 6X = 32 cosh^X-48 cosh^X+18 cosh^X-1 

cosh 7X = 64 cosh^X-112 cosh^X+56 cosh^X-7 cosh X 

Using Table B.l, the first seven polynomials set equal to zero in the 
denominator of equations (B-16) and (B-17) are listed in Table B.3. 

Table B.3 

Expansions of Impedance Function Denominator 



0 = 1- cosh X = 0 (m=l ) 
0=1-2 cosh X (m-2) 
0=1+2 cosh X-4 cosh^X (m=3) 
0=1-4 cosh X-4 cosh^X+8 cosh^X etc. 



0=1+4 cosh X-12 cosh^-8 cosh^X+16 cosh^X 
0=1-6 cosh X-12 cosh^X+32 cosh^X+16 cosh^X-32 cosh^X 
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The equations listed in Table B.3 are the equations which must be solved 
to obtain the poles of L-section models of from one to seven sections. 

In obtaining the equations of Table B.3 the sinh X term is divided out 
since sinh X ^ 0. A more thorough discussion of the constant X is now 
in order. In general, X = + jX^. , so that in equation (B-18), for 

example, 

2| r 
“ L C 

cosh X = cosh X„ cos X. + j sinh X„ sin X. =1 s . (B-19) 

r 1 r 1 c 

Equation (B-19) implies: 

2i r 
w L C 

cosh X^ cos X^ = 1 2 — (B-20) 

sinh X^ sin X^. = 0 (B-21 ) 

But, again, 

2| r 
u L C 

cosh X = cosh X^ cos X^ = 1 - — ^ (B-23) 

and it is seen that cosh X is real, and X is pure imaginary (= jX^.) if 
2i r 

>0. With X = jX^. , X^ = 0 so that finally, 

(/l C 

cosh X = cosh jX^ = cos = 1 2 ^^ • (B-24) 

Equation (B-24) shows that all values of cosh X must be real and lie 
between minus one and plus one. This, in turn implies that all 
solutions to the "equation" in the denominator of (B-17) must also be 
real and lie between minus one and plus one i.e., 

-1 .0 _< [cosh _< 1 .0 . (B-25) 

But the denominator of equation (B-17) is the same as that of 
equation (B-16), so equation (B-25) applies for the lossy line also. 
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Thus, the equations listed in Table B.3 are polynomial equations in the 
real variable cosh X for which there are numerous "root-finding" sub- 
routines in any good computation facility. With the roots of equations 
like those in Table B.3 readily available, equations (B-7) or (B-16) 
can be solved for s or co respectively to give the poles: 

S.c. = - (B-26) 

m = 1 ,2,3.... 



"i.c. ° r Jh n-(cosh 

m = 1 ,2,3 

Where use has been made of the fact that: 
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= — , and = — (B-28) 

s m s m s m 

Equations like those in Table B.3 are seen to be solved easily enough, 
but the problem of how to generate the equations up to fifty or more 
sections (m) remains. The details of the generation problem are left 
to Appendix C. 

B.2 Derivation of L-$ection Model Impedance, Capacitor Input, Shorted 
Termination 

Using the m-section line shown in Figure 3.6, write difference 
equations near the nth section (0 < n < m): 
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Subtracting (B-30) from (B-29): 



Vi ■ * Vi = *-3 ^ (I^ - 

Substitute (B-31) into (B-32): 

d^V 

Vl - " Vl = rr 



dt 



imilar to the development in Appendix B.1.1 assume, 

= lie[v„ e*-*"*] 
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and obtain; 
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cosh X = 1 — ^ 

V = A cosh nX + B sinh X 



To obtain the current, I^, assume = Re[I^ e^^^] in equati 



on 



jo)L I = y - V 
s n n+1 n 

^ 1 

^n “ 3u)i~ cosh(n+l)X + B sinh(n+l)X - A cosh nX - B sinh 
With a short at the zeroth node Vq = 0 implies A = 0, and 

/s 

= B sinh nX 

^ B 

^n ~ jtoL [sinh(n+l)X - sinh nX]. 

Finally, for the driving point impedance at the mth node. 



sinh mX 



« 2 I :> I mi iHA *10 0 

s.c. s sinh(m+l }X-sinh mX “ 1,2,3.... 



(B-32) 



(B-33) 

(B-34) 

(B-29): 

(B-35) 

nX] 

(B-36) 

(B-37) 

(B-38) 

(B-39) 



and using equation (B-33); 
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S.C. i yfc 1>2,3 (B-40) 



The solution for the poles of equation (B-39) is indicated in Appendices 
C.2.2 and D.2. 

B.3 Derivation of L-Section Model Impedance, Open Termination 

Using the m section line shown in Figure 3.3, but with the short 
removed, = 0, and s = jw, the derivation proceeds as in Appendix 
B. 1 . 1 to obtain; 



= A cosh nX + B sinh nX 



I = - 



1 



n 



A[cosh(n-l )X-cosh nX] + B[sinh(n-l)X-sinh nX] 



(B-41) 

(B-42) 



But now, Iq = 0, and equation (B-42) implies. 



B = 



cosh X-1 



sinh X 

Substituting (B-43) into (B-41) and (B-42): 

sinh (n+1 )X-si nh nXj 



Vn = A 



= - 



A 



n jo)L^ 



sinhX 

2 sinh nX-sinh(n+l )X - sinh(n-l)xj 



sinh X 



So that looking into the mth section: 



^0 . c . 



2 s 



sinh(m+l )X-sinh mX ~ [ 

inh mX-sinh(m+l)X-sinh(m-l )xj 



As in Appendix B.1.1 , 



cosh X = 1 - 



2, r 
0) L L 
S S 



from which 



0) 



o.c. Z 



Vfec ■ *Voot 



(B-43) 



(B-44) 

(B-45) 



(B-46) 



(B-18) 



(B-47) 
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The solution for the poles of equation (B-46) is indicated in Appendices 
C.2.3 and D.3. 

B.4 Derivation of ir-Section Model Impedance, Shorted Termination 
Using Figure B.l, write difference equations near the nth 



section: 



dl. 



= L. 



n n+1 ^s dt 

dVn 

n ‘'s dt 

I - I = r ^'*'1 
* n+1 S dt 

Substitute (B-50) into (B-48) to eliminate I*; 

dV 



= L, 






'n+1 



+ I 



n n+1 ‘-S dt I s dt ’ "n+1 

Substitute (B-50) into (B-49) to eliminate I*: 



.]• 



^n ■ ^s dt 



dV„.i dV 

n+1 T _ p n 

■ ^+1 " S dt" • 



(B-48) 

(B-49) 

(B-50) 



(B-51) 



(B-52) 



Write the equivalent to equation (B-51) for the next section toward the 
left: 

d 



''n-l - % ■ 4 3t 



dV -| 

^sdf ^'nj • 



From equation (B-52) write: 



(B-53) 



■n.1 = >n - S ^ • 

Substitute (B-54) into (B-51) to eliminate I : 

n+1 



d^V 



''n - ''n.l 



L.C. 



n+1 



's s dt 
Rewrite equation (B-53): 



+ L 



din h2 

n _ L c d 



s dt 



s-s (''n+l * ''nl- 



V 1 - V 
n-l n 



d^V dl 

L.C. — / + 
s s dt^ s dt 



(B-54) 



(B-55) 



(B-56) 
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. . . / X dl 

Substitute (B-56) into (B-55) to eliminate L-r^ • 

at 

ih 

Vl - * ''n-1 = -r 



dt 



Solve equation (B-57) as in earlier cases to obtain: 



cosh A = 1 - u L C 
s s 



and 



= A cosh nA + B sinh nA . 



To find the current use equation (B-56) assuming solutions 



to obtain: 



I„ = Re[I„ 

V„ » Re[V^ eJ“‘] 



'n “ tVl - ''n 



Expanding equation (B-60): 

^ 1 

^ ^ nA + B cosh nA sinh A]. 

With a short at the zeroth node, Vq = 0 implies A = 0 so that 



and 



= B sinh nA 



z„ = - ju,L, - ^ ' . 

n s cosh nA sinh A 



At the mth node or driving point terminals : 



2 = . juL !».\ , 

m s cosh mA sinh A 

The solution of equation (B-64) for the poles is discussed in 
Appendices C.2.4 and 0.4. The results of the solution are the 
of [cosh ^]y.QQ^ which give the poles of (B-64). Also, using 



(B-57) 

(B-58) 

(B-59) 

(B-60) 

(B-61) 

(B-62) 

(B-63) 

(B-64) 

values 



I 
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equation (B-58), 



0) 



s.c. 



'1 - (cosh X) 



root 



or 



0 ) 



S.c. 



m /2L1 - (cosh 
a V [X 



Equation (B-66) follows from (B-65) since 

L = ^ and C = . 

s m s 2m 

B.5 Derivation of ir-Section Model Impedance Open Termination 
The derivation of the open circuit case is the same as 
the short circuit case through equation (B-61). Now, at the "( 
Iq = 0, and equation (B-61) implies B = 0. Thus 



J = . J_ 

n jtoL 



[A sinh A sinh nX] 



= A cosh nX 



m 



cosh mX 



sinh X sinh mX 



m 



m 



Again, 



0) 



o.c. 



m /2U - Uosh 
A V LC 



(B-65) 

(B-66) 

(B-67) 

that of 
"th node, 

(B-68) 

(B-69) 

(B-70) 

(B-71) 



The solution of equations (B-70), (B-71), and (3-1 ) is discussed in 
Appendices C.2.5 and D.5. 






I 
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APPENDIX C 

SOLUTIONS FOR THE POLES OF THE MODEL IMPEDANCE FUNCTIONS 
C.1 Model Impedance Functions 

The model impedance functions obtained in Appendix B for which 
solutions are sought are shov/n in Table C.l. For the purposes of 
Appendix C, "solutions" refers to those values of cosh X in the 
denominator of the impedance function which result in the poles of the 
function. In Appendix B.1.2 it was shown that the denominator of 
equation (B-17) is a polynomial in the real variable cosh X. In 
addition it was shown that the polynomial had real coefficients and 
that all solutions for [cosh should be such that 

-1.0 £ [cosh <1.0. 

Since the denominator of equation (B-17) is a polynomial in cosh X, 
it follows that the denominator of all the impedance functions is a 
polynomial in cosh X with all the properties of the polynomial in 
equation (B-17). It is the denominator polynomial set equal to zero 
which results in the poles of the respective impedance function. If 
the denominator polynomials [D(cosh X)] can be generated for any 
number of sections, m, then it is a simple matter to call upon root 
routines in the computer library to solve for the roots of the denominator 
equation, D = 0. An inspection of Table C.l and reference to Appendix 
B.1.2 shows that the "generation" of D involves first the generation of 
sinh mX and cosh mX for any m. In Appendix B.1.2 it was shown that 

sinh mX = sinh X[Polynomial in cosh X] (C-1) 



i 

I 




i 
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Model 


Driving Pt. Impedance Function 


L-Section 
Inductor Input 
Shorted Termination 


jo)L^ sinh mX 
^m sinh mX-sinh(m-l )X 

“ " [l-(cosh (B-27) 


L- Sect! on 
Capacitor Input 
Shorted Termination 


jtuL^ sinh mX 

^m sinh(m+l) -sinh m (B-39) 

'W (6-40) 


L- Section 
Open Termination 


jtoL^[sinh(m+l )X-sinh X] 

^m ~ 2 sinh mX-sinh(m+OX-sinh(m-l )X 

*)rootl (6-->') 


TT-Section 

Shorted Termination 


-jtuL^ sinh mX 
^m ~ cosh mX sinh X 

“ " jjh n-(cosh (B-66) 


TT-Section 
Open Termination 


-jioL cosh mX 
^m sinh X sinh mX 

“=?yic 'Voot^ (6-^’> 



Table of Impedance Functions for Which 
Poles are Sought 



Table C.l 




4 



142 



cosh mX = [Polynomial in cosh X] . 


(C-2) 


Denote the polynomial in equation (C-1) by B = B(cosh X) 
equation (C-2) by A = A(cosh X) so that 


and that in 


sinh mX = sinh X B„(cosh X) 

m ' 


(C-3) 


cosh mX = A„(cosh X). 

m ' 


(C-4) 



Then it follows that in all the impedance functions of Table C.l the 
denominator has the form 

D(cosh X) = Linear Combination of A(cosh X) and B(cosh X). 

(C-5) 

A routine must be found which will enable the computer to generate A^^ 

and B„, combine them to form D„ = 0, and then solve for the values of 
m m 

[cosh use in 

“ = f ^>rootJ • 

C.2 Generation of Impedance Function Denominator Polynomial 
In general 

sinh mX = sinh(m-l)X cosh X + cosh(m-l)X sinh X, 

so that in the terminology of Appendix C.l, 

sinh X B„ = sinh X B„ , cosh X + A„ , sinh X, 

m m- I m- 1 

and 

1 cosh X + A , . (C-7) 

m m- 1 m- 1 

Simi larly 

cosh mX = cosh(m-l)X cosh X + sinh(m-l)X sinh X, 



so that 
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or 




(C-8) 



Throughout the remainder of tliis appendix denote 

cosh X by "C" 



and 



sinh X by "S" 

so that equations (C-7) and (C-8) become; 




m = 1,2.... (C-9 ) 



= C A„ , + (C^-1)B„ , m = 1,2.... (C-10) 

m m-1 m-1 ’ ' 



With m = 1,2 comparison of equations (C-9) and (C-10) with Tables 

B.l and B.2 indicates that with Aq = 1 and Bq = 0 equations (C-9) and 
(C-10) form a recursion loop for generating all B^ and A^^. In 
equations (C-9) and (C-10) it is the coefficients of the powers of cosh X 
which are needed to solve for the roots of the "D" polynomial, so that 
in B_ and A it is the coefficients v^hich are of interest. 

In their most general form the polynomials A(m) and B(m) take 
the form (for m = 1,2..,.): 



A(m) = a^^^(m)c"’ + a|^(m)c'^"^ + a^ (m) 



(C-11) 



B(m) = b^(m)c'"-^ + b^_^ (m)C'"-2+ b^(m) . 



(C-12) 



Where 



m = number of sections in the model 



C = cosh X 



A(m) = "A" polynomial when model has m sections 
B(m) = "B" polynomial when model has m sections 
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a^. (m) = coefficient of when model has m sections 

b^. (m) = coefficient of when model has m sections 

The subscripts of the coefficients begin at "1" in order to be compatible 

with the subscripting convention in the Fortran IV language. Note from 

Tables B.l and B.2 that 

B(1 ) = 1 impl ies b-j (1 ) = 1 , 

A(l) = C implies a 2 (l) = 1, (C-13) 

and a-| (1 ) = 0. 

To compute B 2 , for example, use equations (C-9), (C-10), and (C-11): 

B 2 = B(2) = C B(l) + A(l) 

B(2) = C[b^(l)] + [a 2 (l) C + a^(l)] 

B(2) = [a 2 (l) + b^(l)] C + a^(l) 

Where 

[a2(l ) + b^ (1 )] = b2(2) 
and 

a^d) = b^(2) 

In general : 

B(m) = C B(m-1 ) + A(m-1 ) 

= C[b^_^ (m -1 )c"^'^+b^_2('"*’' (m-1 )] 

+ [a_(m-l )c"’"Va , (m-1 a, (m-1)] . 

m iH” I I 

Combining coefficients of like powers: 

B(m) = + * a^(m-l)]c""' 



+ 
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+ [b-| (m-1 ) + a2(ni-l )]C 

+ [a^(m-l)] . (C-14) 

Inspection of equation (C-14) shows that 

b^(ni) = [b^._^(m-l) +a^. (m-1)] 2< i £m (C-15) 

(m) = a^(m-l) . i = 1 (C-16) 

A similar procedure applied for A(m) gives: 

a^(m) = [a^._^(m-l) + b^._ 2 (m-l)] i=m, m+1 (C-17) 

ai*(m) = [a^._i(m-l) + b^._ 2 (m-l )-b^. (m-1 )] i=2,3,. . . ,m-l 

(C-18) 

a 2 (m) = [a^. (m-1) - b 2 (m-l)] i=2 (C-19) 

a-|(m) = [-b^(m-l)] . i=l (C-20) 

Given the initial values for A(l) and B(l), equations (C-15) through 
(C-20) can be implemented in Fortran IV, and the computer made to 
generate and store the coefficients of the "A" and "B" polynomials. 
Having the coefficients of "A" and "B" up to, say, fifty sections 
(m=50), the "D" polynomial can be generated by simple algebra. 

C.2.1 Polynomial for L-Section .'lodel , Inductor Input, Shorted 

For the case of the L-Section, inductor input model with a 

shorted termination, the equation for which the roots are sought is 

(from equation B-17) 

sinh mX - sinh(m-l)X = 0 
sinh X[B(m) - B(m-l)] = 0. 



1 
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The case X = 0 has no meaning for this model and so 

B(m) - B(m-l) = 0, (C-21 ) 

and equation (C-21) is seen to be the relation which gives the poles of 
equation (B-17). Thus 

U(m) = B(m) - B(m-l) = 0. (C-22) 

Equation (C-22) must also be written for implementation into Fortran IV 
in the same way that A(m) and B(m) were implemented. Thus: 

D(m) = d^(m)c'"-'' + d^_^(m)c"’'^ d 2 (m)C + d^(m) 

= B(m) - B(m-1 ) 

= [bj^(ni)c"’"^ + b^_^(m)c"’“^ + b 2 (m)C + b^(m)] 

• I^Vl ^ +....b 2 (m-l) + b^(m-l)]. 

Combining coefficients of like powers: 

b(m) = bjj^(m)c'^~^ + [b^_-|(m) - b|^_-| (m-1 )]C*"'^ + 

t‘>ra-2<'"> - * 

.,.[b 2 (m) - b 2 (m-l)]C + [b-|(m) - b-|(m-l)]. 

So that identifying coefficients of D(m): 

d^. (m) = b^. (m) - b^. (m-1) i = l....m-l (C-23) 

d„(m) = b„(m) . (C-24) 

m m 

Thus the computer can be made to generate the coefficients of the "D" 
polynomial. The form of the coefficients given in equations (C-24) and 
(C-25) are explicitly suited for several "root" subroutines available 
in the IBM- 360 Scientific Subroutine Package (SSP) at the Computation 
Center of M.I.T. Equations (C-15) through (C-20) and (C-23) and (C-24) 



i 
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were implemented in Fortran and solved with the IBII-360 as shown in 
Appendix U.l. In Appendix C.l the roots of L) are put into equation 
(B-26) to obtain the poles of a 100 mile long typical 765 kV line. 

The poles of the model thus obtained are compared to the poles of the 
actual line to obtain the pole error plots of Chapter 3. 

C.2.2 Polynomial for L-$ection llodel , Capacitor Input, Shorted 

For the case of the L-section, capacitor input model with a 
shorted termination, the equation for which the roots are sought is 
(from equations B-39): 

sinh(m+l)X - sinh mX = 0 
sinh X[B(m+l) - B(m)] = 0. 

The case X = 0 has no meaning for this model and so 

D(m) = B(m+1) - B(m) = 0. 

Following a procedure identical to that in Appendix C.2.1, the 
coefficients of the "D" polynomial are found to be 

cim(m) = bj^(m) (C-25) 

d^(m) = [b^. (m)-b^ (m-1 )]. i = l ,2, . . . .m-1 (C-26) 

iJote that equations (C-25) and (C-26) appear identical to equations 
(C-24) and (C-23) respectively, but that in the case of equations 
(C-25) and (C-26) the coefficients pertain to a higher power of the 
variable C. Equations (C-15) through (C-20) and equations (C-25) 
and (C-26) were implemented in Fortran and solved with the computer 
as shown in Appendix 0.2. 
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C.2.3 Polynomial for L-Section Model, Open Termination 

For the case of the open circuited L-section model, the equations 
for which the roots are sought is (from equation B-46); 

2 sinh mA - sinh(m+l)A - sinh(m-l)A = 0 
sinh mA(l - cosh A) = 0 
sinh A B(m)[l-C] = 0. 

Again the case A = 0 adds no solutions, so that division by sinh A can 
be performed assuming A 0. Finally then 

D(m) = B(m)(l-C) = 0 

D(m) = [b^(m)c"’’'' + b^_^(m)c'"’" +....b 2 (m) + b^(m)](l-C). 

Collecting coefficients of like powers of C: 

D(m) = - b^(m)c"' + [b^(m) - (m)]c"’'^ + 

+[b 2 (m) - b^(m)]C + b^(m) . 

The coefficients, d^. , are then seen to be 

d^(m)=b^(m) (C-27) 

dj(m) = [b^.(m) - b^._-|(m)] i=2,3, m (C-28) 

dm+i(m) = - b|^(m) when i=m+l (C-29) 

Equations (C-15) through (C-20) and (C-27) through (C-29) were 
implemented in Fortran and solved with the computer as shown in 
Appendix D.3. 

C.2.3 Polynomial for n-Section Model, Shorted 

For the case of the shorted Tr-section model, the equation for 
which roots are sought is (from equation B-64): 



cosh mA sinh A = 0. 
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The case A = 0 has no meaning for this model so that 

D(m) = cosh mA, 

0(m) = A(m), 

d^. (m) = a^(m) i = l,2,....m+l 

This case is particularly simple and is shown implemented in Appendix D.4. 
C.2.4 Polynomial for Tr-Section liodel , Open 

For the case of the open ir-section model, the equation for which 
roots are sought is (from equation B-70): 

sinh A sinh mA = 0 
sinh^A B(m) = 0 
(1-cosh^) B(m) = 0 

D(m) = (1-C^) B(m) = 0 . (C-30) 

Expanding equation (C-30) and collecting like powers of C: 

0(m) = b^(m)c”*' + 6 „.,(in)c"’ t [b^.2(m) - b^(ai)]c"''’ +... 

[b2(m) - b^(m)]C^ * [b,(m) - bj(m)]C^ - b2(m)C - b,(ra). 

The coefficients d^- are then identified as 



d.(m) = - b. (m) if i 


= 1,2 


(C-31) 


d^.(m) = [b^._2(ni)-b^. (m)] 


i=3,4,. . .m 


(C-32) 


d. (m) = 


i =m+l 


(C-33) 


di(m) = hjm) 


i =m+2 


(C-34) 



Equations (C-15) through (C-20) and (C-31) through (C-34) were 
implemented in Fortran and solved with the computer as shown in 
Appendix 0.5. 
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APPENDIX D.l 
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